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During  this  report  period,  the  followin'-  people  have  been  involved  in  the 


exper iment : 


distances . 


Since  room  temperature  gradiometers  are  ultimately  limited  by  room  tempera¬ 
ture  Brownian  motion,  we  have  been  experimenting  with  a  laboratory  cryogenic 
gradiometer  which  could  lead  to  the  development  of  a  more  sensitive  moving  base¬ 
line  gradioraeter  for  field  use.  Such  gradiometers  could  have  applications  to 
problems  in  navigation  and  modeling  of  the  earth's  gravitational  field. 

Recently  the  inverse  square  law  of  gravity  at  laboratory  distances  has 
become  a  subject  of  great  interest  to  physicists.  This  has  come  about  for  two 
reasons.  First,  Long^  and  more  recently,  Tuck^  have  found  an  apparent  viola¬ 
tion  to  the  inverse  squire  lav  of  ;\i  'vlty  at  laboratory  dist  ices.  If  this  were 
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William  Fairbank  -  Co-Principal  Investigator 

t  | 

Daniel  DeBra  -  Co-Principal  Investigator  fljg SWJ 

Michael  McAshan  -  Senior  Research  Associate 

Evan  Mapoles  -  Graduate  Student  and  then  Research  Associate.  His  Ph.D. 
thesis  entitled  "Development  of  a  Superconducting  Gravity  Gradiometer  for  Test 
of  the  Inverse  Square  Law,"  is  included  as  Appendix  1.  His  publication, 

"A  Superconducting  Gravity  Gradiometer"  is  Included  as  Appendix  2. 

RESEARCH  OBJECTIVES  ’  < 

r  -* 
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This  research  has  had  two  objectives:  (1)  to  study  the  application  of  •' 

cryogenic  technology  to  the  problem  of  gravity  gradient  measurements,  and 
(2)  to  use  such  a  gradiometer  to  make  more  accurate  meaurements  of  the  inverse 
square  law  of  gravity  to  test  the  recent  experimental  and  theoretical  sugges- 
tlons  that  the  inverse  square  law  of  gravity  might  be  violated  at  laboratory  .*‘/.va 
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If  this  were 


found  to  be  correct,  ft  would  be  one  of  the  most  important  experiments  in  the 


history  of  physics.  Secondly,  theorists  have  lugun  to  speculate  on  the  possible 

existence  of  new  fundamental  particles  which  could  affect  the  inverse  square  law 

of  gravity  at  laboratory  distances.  For  example,  axions  with  mass  10  ^  eV  have 

been  postulated  to  explain  the  dark  matter  of  our  galaxy  and  the  reason  for  the 

lack  of  parity  in  time-reversal  violation  and  strong  interactions.  Such  particles 

(3) 

could  contribute  to  a  violation  of  the  inverse  square  law.  By  making  use  of 

a  superconducting  gradiometer  it  appears  possible  to  check  the  inverse  square 
law  at  distances  of  the  order  of  one  meter  to  an  accuracy  of  2  parts  in  10^. 

Operation  of  a  gravity  gradiometer  at  4  K  offers  a  reduction  of  thermal 
noise  power  by  a  factor  of  75  over  that  in  room  temperature  environments,  and 
this  is  the  principle  motivation  for  the  development  of  cryogenic  instruments. 
However,  there  are  many  other  potential  advantages  which  help  to  offset  the 
additional  effort  required  to  maintain  the  instrument  at  4  K.  Below  10  K 
several  metals  lose  all  low  frequency  electrical  resistance  (superconductivity) 
and  become  nearly  perfect  magnetic  shields.  These  materials  are  used  In  our 
gravity  grauiometers  to  make  magnetic  readouts  which  are  free  of  thermal  noise. 
Although  careful  construction  is  required  to  allow  for  the  thermal  contraction 
of  the  instrument  as  it  is  cooled  from  room  temperature  to  4K,  once  the  device 
is  cold  the  coefficient  of  thermal  expansion  is  far  below  the  room  temperature 
value.  The  liquid  helium  bath  commonly  used  to  maintain  low  temperature  also 
provides  a  stable  thermal  environment.  In  addition,  low  frequency  amplifiers 
of  unequaled  sensitivity  based  on  the  Josephson  effect  (SQUIDs)  are  commercially 


ava l Table. 


SUMMARY  OF  RESEARCH  EFFORT 

During  the  early  part  of  our  work  we  constructed  and  tested  two  prototype 
gradiometers.  These  devices  both  incorporated  superconducting  coils  and  proof 
masses,  were  built  with  similar  techniques,  and  both  measured  a  single  diagonal 
component  of  the  gravity  gradient  tensor.  They  differed  in  the  method  used  to 
obtain  common  mode  balancing.  One  of  our  prototypes  used  four  superconducting 
coils,  one  on  each  side  of  two  superconducting  diaphragms.  Each  pair  of  coils 
and  diaphragm  produced  a  current  in  response  to  an  applied  acceleration,  and 
these  two  currents  were  subtracted  at  the  input  of  a  SQUID.  We  referred  to 
this  gradioneter  as  "current  differencing".  Our  second  prototype  used  a  single 
coil  located  between  two  superconducting  proof  masses  so  that  its  inductance 
was  modulated  only  by  relative  notion  of  these  two  proof  masses.  We  called 
tills  gradiometer  "displacement  differencing".  Because  the  four  coils  in  the 
current  differencing  gradiometers  are  not  perfectly  matched,  the  overall 
balancing  problem  in  this  gradiometer  is  more  complicated.  Consequently,  we 
choose  the  displacement  differencing  design  to  an  improved  version. 

We  have  completed  construction  of  this  improved  gradiometer  and  operated 
it  at  4.2  K.  with  sensitivity  approaching  1E/-7  Hz  at  frequencies  between  2  and 
10  II.’..  The  operation  of  this  gradiometer  is  described  in  Appendices  1  and  2. 
in  order  to  achieve  this  level  of  sensitivity  at  lower  frequencies,  a  new  cryo¬ 
stat  probe  was  built  to  provide  a  1.3  K  environment  for  the  gradiometer.  The 
purpose  of  this  was  to  reduce  the  thermal  sensitivity  of  the  Instrument.  In 
addition,  the  new  probe  was  designed  to  provide  a  high  degree  of  seismic  and 
magnetic  isolation.  Thin  probe  Is  being  used  to  gather  data  from  the  gradio¬ 
meter  in  the  Inverse  square  It..’  experiment. 
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DESCRIPTION  OF  THE  GRAD IOMETER 


Our  primary  goal  in  the  design  and  fabrication  of  our  gravity  gradiometer 

was  to  minimize  the  instrument's  sensitivity  to  linear  accelerations  both  along 

the  proof  mass  direction  and  perpendicular  to  it.  Linear  accelerations  in  the 

laboratory  exert  forces  on  the  proof  masses  which  are  very  large  compared  to 

the  gravitational  forces  we  would  like  to  measure,  so  careful  balancing  of  the 

gradiometer  is  needed  to  eliminate  them  as  a  noise  source, 

In  addition,  we  decided  to  build  an  on-line  gradiometer  which  measures  one 

of  the  diagonal  components  of  F.  It  is  highly  desirable  to  be  able  to  measure 

..,2 

all  three  components  of  v  J1  so  wc  chose  a  method  for  suspending  the  proof  masses 
which  allows  us  to  operate  the  gradiometer  with  its  sensitive  axis  vertical  as 
well  as  horizontal. 

In  order  to  make  the  balancing  problem  manageable  we  chose  a  mechanical 
design  which  enabled  us  to  use  only  one  coil  in  the  readout  system  described 
later.  Experiments  have  been  done  with  readouts  using  several  coils,  but  in 
these  gradiometers  careful  matching  of  the  coils  as  well  as  the  proof  masses 
is  necessary  to  reject  linear  accelerations. 

Basically  the  gradiometer  consists  of  two  superconducting  proof  masses 
adjacent  to  each  other.  A  single  niobium  coil  Is  mounted  In  the  gap  between 
the  proof  masses  so  that  the  coil  inductance  is  proportional  to  their  spacing. 

This  readmit  system  is  analyzed  in  detail  below.  Since  the  sensitive  axes  of 
the  proof  masses  are  parallel,  the  gradiometer  detects  changes  in  one  of  the 
diagonal  components  of  the  gravity  gradient  tensor,  T 

j  j 

Each  of  the  two  proof  masses  is  supported  with  two  mechanical  springs. 

I  he  mechanical  springs  are  folded  c  ini'  levers  cut  into  ci’-'-ulnr  di^'-.s  of 
niobium.  A  foM-  !  oautil'ver  viih  tvo-h'ld  symmetry  is  shown  in  figure  1. 


In  this  case  the  thickness  of  the  plate  <  s  reduced  in  the  two  rectangular 
areas  hatched  in  the  figure  in  order  to  restrict  the  bending  to  those  areas. 

Slots  are  then  cut  in  the  indicated  locations.  These  slots  allow  the  center 
to  move  perpendicular  to  the  plane  of  the  plates  shown  in  the  lower  part 
of  the  figure.  Since  the  center  of  the  spring  is  free  to  move  in  the  plane 
of  the  plate,  the  total  length  of  the  spring  remains  constant  and  the  restor¬ 
ing  force  is  provided  by  pure  bending  without  stretching  of  the  spring  material. 

The  proof  mass  subassembly  consists  of  two  springs,  two  cover  pieces,  an 
annular  mass,  and  a  ring  used  to  hold  the  edges  of  the  springs  apart.  The  cover 
pieces  are  circular  disks  of  niobium  8.51  cm  in  diameter  which  are  flat  on  one 
side  and  have  a  threaded  stud  on  the  opposite  side.  The  flat  sides  are  designed 
to  modulate  the  inductance  of  the  readout  coil  or  tuning  coils  described  later. 

The  flats  were  polished  to  ensure  flatness  and  Isotropy  so  that  close  spacing 
to  the  colls  could  be  obtained,  and  so  that  the  inductance  of  the  coils  is  modu¬ 
lated  only  by  linear  motion  of  the  surfaces.  The  threaded  studs  screw  into  an 
annular  disk  of  niobium  which  provides  the  bulk  of  the  mass  for  each  proof  mass. 

A  spring  is  captured  between  each  cover  piece  and  proof  mass.  An  exploded  view 
of  one  of  the  two  proof  mass  subassemblies  is  shown  in  Figure  2.  The  ring  which 
is  used  to  maintain  the  spacing  of  the  edges  of  the  springs  is  thicker  than  the 
proof  mass  so  that  each  spring  is  biased  away  from  its  equilibrium  position  and 
never  passes  through  its  zero  point .  This  eliminates  certain  forms  of  nonlinear 
behavior.  Since  the  springs  have  opposite  biases,  the  first  order  nonlinear  terms 
also  cancel.  1  he  ring  was  grounded  to  its  final  thickness  to  keep  its  two  sides 


parallel  to  one  part  in  1U  . 


Once  the  two  proof  mass  subassemblies  were  constructed,  they  were  stacked 
on  top  of  each  other  with  a  ring  holding  them  apart.  The  thickness  or  the  ring 
was  chosen  so  that  the  readout  coil  is  about  0.013  cm  from  the  opposite  proof 
mass.  This  structure  was  mounted  on  an  aluminum  fixture  which  maintained  the 
angular  orientation  of  the  parts.  A  niobium  housing  tube  9.3  cm  long  with  a 
10.5  cm  inside  diameter  was  placed  over  the  parts  (Figure  3).  A  more  detailed 
description  of  this  grad  Lumet  or  is  given  in  Appendix  1,  including  details  of 
caper  imental  measurements  on  the  ".radiometer. 

The  assembled  grniliomoler  was  first  tested  at  4.2  K  and  has  performed  as 
expected  i  u  tin'  frequency  range  from  2  !b:  to  1.0  Ilx  (see  Appendix  1).  It  has 
since  been  Incorporated  in  a  new  cryostat  probe  capable  of  cooling  the  instru¬ 
ment  to  1.3  K  and  is  being  used  on  the  inverse  square  law  experiment. 

K.Vi'RSi:  FOi’ARC  i. AW  KXTF.K  I  MiwiT 

The  theory  and  background  for  this  experiment  is  given  in  Chapters  4  and  5 
of  Appendix  1.  We  have  completed  construction  of  all  major  components  of  the 
inverse  square  law  experiment  (Figure  4)  and  have  commenced  work  on  the  experi¬ 
ment.  Those  components  include  a  1.3  K  cryogenic  probe  and  a  pedestal  which  is 
usd  to  mount  Lhe,  cryostat  in  a  ton  foor  hole  in  the  floor.  A  pumping  line  is 
incorporated  into  the  pedestal  which  allows  us  to  continuously  pump  on  the 
cryostat,  to  maintain  the  temperature,  of  Lhe  gradiometer  at  1.3  K  while  the 
cylindrical  source  is  moved  up  and  down  over  the  cryostats.  The  final,  major 
pins,'  of  the  apparatus  is  Lhe  .system  for  moving  lhe  cyLinder  up  and  down. 
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FIGURE  CAFTIONS 


Assembled  gradlometer  in  mounting  fixture.  Mounting  fixture 
includes  flex  hinges  on  both  sides  of  the  gradlometer  so 
that  the  instrument  may  be  tilted  slightly  by  PZT  bender 
plates  (upper  center). 

Major  components  of  the  inverse  square  law  experiment.  The 
cylindrical  test  mass  is  shown  in  the  upper  and  central  portion 
of  the  picture  and  the  dewar  on  its  pedestal  is  barely  visible 
at  the  bottom. 

Upper  and  mid-portion  of  the  new  cryostat  probe.  Cold  plate 
and  the  upper  part  of  the  1.3  K  radiation  shield  are  at  the 
very  bottom.  The  isolation  stack  (which  would  be  at  the  very 
top)  is  not  included  at  this  level  of  assembly. 

Lower  detail  of  the  1.3  K  cryostat  probe  including  the  gradio- 
meter  and  its  fixture  in  the  lower  portion  of  the  picture.  The 
plate  in  the  central  portion  is  a  1.3  K  continuous  cold  plate 
refrigerator.  During  operation  the  gradiometer  is  completely 
enclosed  by  copper  radiation  shield  attached  to  the  cold  plate. 
In  turn  the  entire  lower  portion  od  the  probe  is  isolated  from 
the  4  K  helium  bath  by  a  vacuum  can  sealed  onto  the  top  plate 
seen  in  the  upper  part  of  the  picture. 

Mechanical  isolation  stack.  This  assembly  sits  on  top  of  the 
cryostat  and  provides  the  sole  mechanical  support  for  the  gradio 
meter.  Durim;  operation,  the  stack  is  enclosed  in  a  vacuum  can. 


In  order  to  maintain  the  gradiometer  at  1.3  k,  the  new  probe  has  a  copper 
plate  inside  of  the  vacuum  can  which  serves  as  a  refrigerator  (Figure  5).  The 
probe  vacuum  can  is  immersed  in  a  liquid  helium  bath  at  4.2  K.  The  copper 
plate  inside  the  can  is  thermally  isolated  from  the  can  and  contains  a  small 
reservoir  which  is  connected  to  the  4.2  K  bath  through  a  flow  impedance  (Figure  6). 
This  reservoir  is  continuously  pumped  to  maintain  the  plate  at  1.3  K.  The  flow 
impedance  is  adjusted  at  room  temperature  to  adjust  the  temperature  of  the  plate. 
The  gradiometer  is  thermally  grounded  to  this  plate  to  operate  it  at  1.3  K.  The 
cold  plate  has  been  tested  successfully.  It  has  the  major  advantage  that  liquid 
helium  can  be  added  to  the  bath  in  the  usual  way  without  disturbing  the  1.3  K 
space.  This  means  that  the  1.3  K  space  can  be  maintained  for  long  periods. 

A  second  major  improvement  in  the  new  probe  is  the  addition  of  a  mechanical 
transmission  line  which  will  isolate  the  gradiometer  from  vibrations  at  frequen¬ 
cies  above  10  Hz  (Figure  7).  This  line  consists  of  four  brass  disks  which  hang 
from  each  other  with  small  latex  springs.  The  soft  springs  cut  off  vibrations 
above  a  frequency  of  about  10  Hz.  The  gradiometer  hangs  from  the  bottom  of  this 
stack  by  a  stainless  steel  tube  which  allows  the  gradiometer  to  be  at  1,3  K 
while  the  stack  is  at  room  temperature.  The  entire  assembly  is  in  vacuum  to 
isolate  from  acoustic  noise. 

In  order  to  measure  the  gradient  at  the  center  of  the  cylindrical  gravita¬ 
tional  .source,  the  dewar  containing  the  probe  and  the  gradiometer  has  been  mounted 
on  a  pedestal  which  will  fit  inside  of  the  cylinder.  The  pedestal  is  mounted 
at  the  bottom  of  a  ten  foot  deep  hole  which  the  cylinder  is  lowered  into  during 
the  measurement .  ibis  leave:;  the  dewar  near  floor  level  to  facilitate  adjust¬ 
ments  .  A  pumping  line  i  s  built  in  so  that  the  cold  plate  can  be  pumped  contin¬ 
uously  as  t  lie  cyl  inder  goes  up  and  ilo.-n.  The  7,160  pound  cvl  imlcr  rolls  along 
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lo  stool.  I-beam  up  and  down  the  wall  of  the  end  station  in  which  the  experiment 
takes  place.  A  hoist  has  been  installed  to  raise  and  lower  it.  Considerable 
care  has  been  taken  to  ensure  smooth  motion  of  this  large  mass,  and  to  minimize 
the  vibration  which  the  gradiometer  must  tolerate  during  the  experiment. 

During  the  past  year  the  completed  gradiometer  and  1.3  K  cryostat  assembly 
has  been  subjected  to  initial  tests.  During  these  tests  gradiometer  output  power 
spectra  were  obtained  with  the  gradiometer  in  unbalanced  operation.  (Balanced 
operation  was  not  possible  due  to  inability  to  store  sufficient  current  in  the 
balancing  circuit.  This  problem  is  presumed  to  be  due  to  ohmic  heating  in  the 
current  supply  leads  whirl)  have  been  subsequently  improved.)  The  output  power 
spectrum  obtained  at  1.3  K  was  observed  to  have  two  distinct  regions.  Above  about 
0.2  Hz  the  spectrum  was  dominated  by  externally  driven  mechanical  motion  of  the 
gradiometer,  most  notably  the  pendulum  modes  of  the  suspension  system.  Previous 
experience  with  the  gradiometer  has  indicated  that  this  response  is  reduced  by 
a  factor  of  -1000  when  the  gradiometer  is  properly  balanced.  Below  ~0.1  Hz 
the  power  spectrum  was  dominated  by  low-frequency-divergent  behavior  arising 
from  the  temperature  coefficient  of  the  gradiometer  together  with  temperature 
fluctuations  of  the  device.  In  the  frequency  range  of  interest  for  the  inverse 
square  law  experiment  (~0.15  I!z)  the  noise  power  is  estimated  to  30-40  E/Viiz 
where  a  calibration  factor  obtained  from  operation  at  4  K  is  used.  Proper 
balancing  of  the  gradiometer  should  drop  this  to  a  SQUID  noise  limited  value 
of  1  E/*Hz.  Further  reduction  of  this  value  by  an  order  of  magnitude  seems 
possible  by  the  use  of  a  dc  SQUID  which  is  now  commercially  available. 


Although  the  maximum  frequency  at  which  temperature  fluctuations  dominate 
the  noise  power  spectrum  has  been  shifted  downward  from  1  Hz  at  4  K  to  0.1  Hz 
at  1.3  K,  the  temperature  coefficient  at  the  lower  temperature  is  higher  than 
expected.  In  addition,  temperature  fluctuations  at  the  gradiometer  are  on  the 
order  of  10  mK,  much  larger  than  what  is  possible  to  obtain.  We  have  tenta¬ 
tively  identified  a  free-floating  radiation  baffle  located  on  the  gradiometer 
suspension  rod  between  the  cold  plate  and  the  top  of  the  vacuum  can  as  a 
possible  source  of  this  problem.  We  are  in  the  process  of  mounting  a  heater 
and  thermometer  on  this  baffle  in  order  to  be  able  to  check  this  hypothesis 
and  will  be  able  to  do  so  as  soon  as  a  leaking  low-temperature  feedthrough 
seal  is  replaced. 

FUTURE  DEVELOPMENTS 

In  the  following  year  this  researcli  will  be  continued  witii  the  aid  of  a 
fourth-year  graduate  student,  Joel  Parke,  who  will  be  visiting  with  Professor 
Ho  Juii},  ”aik  from  the  University  of  Maryland.  Our  goal  is  to  check  the  inverse 
square  law  to  the  highest  possible  accuracy  and  to  understand  and  control  the 
temperature  coefficient  which  remains  at  1.3°  K. 


-*  J1 
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ABSTRACT 

A  gravity  gradlometer  has  been  developed  which 
uses  two  superconducting  proof  Basses  to  measure  a 
single  in  line  component  of  the  gravity  gradient 
tensor .  A  niobium  coll  positioned  between  the  proof 
masses  carries  a  persistent  current  which  is  modu¬ 
lated  by  the  relative  aotlon  of  the  masses.  This 
current  signal  Is  detected  by  a  SQUID.  We  expect 
to  achieve  a  sensitivity  of  10~9/b~2  Hz"**.  Initial 
experiments  indicate  that  the  detector  is  working 
properly.  An  experiment  to  test  the  Inverse  square 
law  of  gravity  at  laboratory  distances  Is  planned. 

I.  INTRODUCTION 

During  the  past  ten  years  a  cryogenic  transducer 
was  developed  at  Stanford  to  convert  a  mechanical 
signal  in  an  aluminum  bar  to  an  electrical  signal  for 
amplification  (Palk,  1976).  This  transducer  is  cur¬ 
rently  used  in  the  Stanford  gravitational  radiation 
experiment.  The  transducer  is  basically  an  accelero¬ 
meter  whose  sensitivity  is  optimized  for  detection 
of  the  fundamental  longitudinal  mode  of  the  bar. 

Since  1975  we  have  been  investigating  the  poss¬ 
ibility  of  using  this  technology  for  gravity  gradient 
measurement.  Potential  advantages  Include  reduced 
thermal  noise  due  to  low  temperature  and  higher 
mechanical  Q’s  of  materials,  increased  mechanical 
stability,  and  the  use  of  SQUID's  as  very  low  noise 
amplifiers.  Two  prototype  gradiometers  using  diff¬ 
erent  methods  of  common  mode  rejection  were  assembled 
with  components  of  gravity  wave  transducers  (Paik, 
19,'d).  The  operating  principles  of  both  svstems 
* 
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were  confirmed  experimentally,  and  one  of  the  methods 
of  common  mode  rejection  was  chosen  for  an  advanced 
version. 

In  the  past  eighteen  months  we  have  designed 
and  built  this  new  gradloneter.  It  Is  a  mechanically 
improved  system  designed  for  higher  conmon  mode 
rejection.  Testing  began  in  April,  1980.  The  mech¬ 
anical  and  electrical  design  of  this  gradlometer,  as 
well  as  some  preliminary  test  results  are  discussed 
here. 

II.  MECHANICAL  CONSTRUCTION 

The  gradlometer  consists  of  two  proof  masses 
held  adjacent  to  each  other  with  their  sensitive 
axes  along  a  line.  The  device  measures  a  single 
diagonal  component  of  the  gravity  gradient  tensor. 
These  two  proof  masses  are  carefully  matched  and 
aligned  to  minimize  sensitivity  to  common  accelera¬ 
tion.  The  sensor  package  is  a  cylinder  9.5  cm  long 
and  11.5  cm  in  diameter.  Each  proof  mass  weighs 
1.07  kg. 

In  order  to  achieve  careful  alignment  of  the 
proof  mass  sensitive  axes,  the  components  of  the 
gradlometer  are  designed  to  stack  up  on  each  other. 
The  faces  of  the  parts  are  made  as  parallel  as  poss¬ 
ible.  This  stack  slides  into  a  niobium  housing  tube, 
where  it  is  locked  together  with  two  end  rings  which 
thread  into  the  housing  tube.  With  this  method  we 
expect  the  axes  to  be  aligned  to  a  few  tens  of  micro- 
radians.  Accelerations  perpendicular  to  the  gradie- 
raeter  axis  cause  an  error  which  is  proportional  to 
this  angular  misalignment.  Consequently  we  expect 
to  reject  these  cross  accelerations  to  a  few  parts 
in  105. 

Each  proof  mass  consists  of  an  annular  disc 
with  a  threaded  central  hole,  two  flat  springs  and 
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two  ground  plane  pieces  which  screw  Into  either  aide 
of  the  annular  dlac.  The  center  of  each  aprlng  Is 
locked  between  the  center  dlac  and  a  ground  plane. 

The  entire  assembly  la  niobium  ao  that  each  proof 
mass  has  two  superconducting  surfaces  which  can  bo 
used  to  modulate  readout  colls.  The  rltns  of  the  two 
springs  are  held  apart  with  a  niobium  ring  which  Is 
fixed  to  the  housing  tube. 

The  flat  springs  are  made  from  solid  discs  of 
niobium  .318  cm  thick  and  10.48  cm  in  diameter.  A 
rectangular  depression  Is  cut  into  the  disc  in  three 
places  120°  apart  so  Chat  it  is  .071  cm  thick  in  an 
area  3.8  cm  by  2.2  cm.  Thi6  area  was  cut  out  with 
an  electron  discharge  machine.  After  the  depressions 
were  cut,  a  aeries  of  very  narrow  (.024  cm)  6lots 
were  cut  with  a  wire  EDM  machine.  With  these  slots 
each  rectangular  depression  becomes  a  folded  canti¬ 
lever  spring.  The  center  of  the  disc  is  thus  weakly 
held  with  respect  to  the  rim.  The  center  is  attached 
to  the  proof  mass  and  the  rim  Is  fixed  to  the  case. 

The  two  proof  mas6  subassemblies  are  held  apart 
Inside  the  housing  tube  with  a  third  niobium  ring. 

In  t lie  gap  between  the  two  adjacent  superconducting 
ground  planes  a  superconducting  pickup  coll  16  moun¬ 
ted.  This  coil  is  used  to  detect  dltferentlal  motion 
of  the  proof  masses.  Two  additional  coils  are 
mounted  facing  the  outer  surfaces  of  the  proof  masses 
for  detection  of  common  motion.  These  also  serve  as 
magnetic  springs  which  can  be  used  to  tune  out  any 
sensitivity  to  common  acceleration  along  the  pradio- 
meter  axis,  figure  1  shows  a  cross  section  of  the 
grad loroeter . 

Except  for  the  thin  dielectric  forms  on  which 
the  coils  are  mounted  the  entire  assembly  is  pure 
niehiur.  Since  niobium  undergoes  a  strain  of  1.43 
x  !■.>  '  when  cooled  from  room  temperature  to  4.2°  k, 
all  parts  must  be  of  the  same  material  to  prevent 
differential  contraction  from  altering  the  geometry. 

In  order  to  calibrate  the  system,  the  gradlo- 
meter  package  is  mounted  on  a  fixture  which  allows 
it  to  be  tilted  by  applying  currents  to  two  super- 
conduit mp  colls  mounted  in  speaker  magnets.  The 
b' radiometer  is  held  bv  two  cross-flex  pivots  perpon- 
Jiiuljr  to  its  sensitive  axes.  The  speaker  coils 
•H»r  1  v  ,i  torque  about  this  axes.  If  a  tilt 
0  *  0  S  i  nil  t 


Is  applied  to  the  grad lometer ,  on  acceleration 


gradient 


Aa  -  u>  20  2  ,  , 

o  o  cos2u  t 


results.  By  measuring  the  amplitude  of  the  second 
harmonic  of  the  drive  signal  the  instrument  scale 
factor  can  be  determined. 


Ci  at  tDDibir 


COIL  FORMS 
Figure  1.  Gradiometer  Cross  Section 

III.  ELECTRICAL  CIRCUITS 
The  gradiometer  contains  two  electrical  circuit 
one  tor  detection  to  differential  proof  mass  motion 
and  one  for  common  motion.  Each  is  readout  with  a 
SQUID. 

The  differential  circuit  is  shown  in  Figure  2. 
Thr  coll  mounted  bet veto  the  proof  masses  Is  a  singl 
layer  of  .0089  cm  niobium  wire  wound  in  a  spiral  or> 
a  .254  cn.  wafer  of  machinable  glass  ceramic.  This 
wafer  is  glued  to  one  of  the  proof  masses  such  that 
the  surface  of  the  wire  is  ver ,  close  (.013  cm)  to 
the  other  proof  mass.  Since  mignetii  field  1 o\,  lu- 
ded  f  rom  the  -uiperrondix  t  i  np  proof  m.is^.,  the  inju«  • 
tancc  of  tin-  pii  kuo  <oi!  i*  *ont  iolle-1  bv  tin  ope¬ 
ration  of  tnc  prooi  mas.*--  wh  ■  ,  ihamev  in  »  e -.-i  >n«-e 


Co  gradients.  A  persistent  current  1  It  stored  In 

a  superconducting  loop  containing  the  pickup  coil 

inductance  L  and  is  coupled  to  a  SQUID  via  a  auper- 
P 

conducting  tranaioraer.  Since  the  flux  In  this  loop 
la  fixed,  as  the  proof  aass  separation  changes  a 
current  i  is  generated  In  the  SQUID  Input  coll.  _ 


I 


Figure  2.  Differential  Mode  Detection 
Circuit 


For  an  Ideal  transformer  the  detection  sensiti¬ 
vity  can  be  calculated  from 


where  1  Is  the  signal  current  produced  In  the  SQUID 
by  an  acceleration  gradient  G,  b  la  the  proof  Bass 
separation,  d  is  the  spacing  of  the  pickup  coll  from 
the  other  proof  mass,  iv  Is  the  resonant  frequency 
of  the  differential  mode,  and  Lp  and  Ls  are  the 
pickup  cml  and  SQUID  Input  coil  Inductances.  In 
order  tc  optimize  the  sensitivity  the  gap  d  is  made 
as  small  as  practical,  and  the  resonant  frequency 
is  minimized.  Because  it  ls  difficult  to  make 
stable  low  frequency  suspensions,  we  chose  60  Hz 
foi  u.  This  allowed  us  to  use  a  practical  thickness 
for  our  folded  cantilever  springs.  With  I  *  8A , 
tiio  S'H  IU  noise  gives  a  gradient  noise  of 
1  K/w  H>  (It  *  1C)"9  S~;  ) . 

This  type  of  readout  is  sensitive  to  temperature 
Decause  the  penetration  depth  of  the  super¬ 
conducting  proof  masses  varies  with  temperature. 

As  flux  moves  into  the  superconductor  the  coil 
inductance  varies  as  If  the  surface  were  moving 
phvsica 1 1 • .  The  dependence  of  penetration  depth 
*  on  temperature  is  given  by 


(Tinkha®,  1975)  where  X  •  440  A  arid  T  »  9.2°  K 
o  C 

for  niobium.  Taking  the  derivative  with  reapect  to 
temperature  we  find  that  at  4.2°  K 


dX 

dT 


9.7  A/°K 


The  differential  displacement  caused  by  s  1  E  grad¬ 
ient  ls  about  2  x  10  ^  a.  Thus  we  Bust  control  the 
temperature  to  10  6  °K  to  stabilize  the  gap.  This 
requirement  can  be  relaxed  somewhat  by  operating 
the  gradioioeter  at  a  lower  temperature. 

The  circuit  shown  In  figure  3  is  used  both  to 
detect  common  motion  of  the  proof  masses  and  to 
fine  tune  the  stiffness  of  the  proof  mass  suspen¬ 
sions.  The  principle  of  detection  ls  the  same  as 
that  uaed  In  the  differential  circuit.  Each  proof 
mass  modulates  a  persistent  current  in  a  flat  coll 
facing  its  outer  surface.  The  error  signals  genera¬ 
ted  by  this  motion  are  gummed  in  a  toroidal  trans¬ 
former  and  coupled  to  the  input  of  a  SQUID. 

The  field  generated  by  the  persistent  currents 
provide  a  restoring  force  between  the  proof  mass  and 
the  case  when  the  mass  ls  displaced  from  equilibrium. 
Hence,  each  pickup  coil  acts  as  a  magnetic  spring 
In  parallel  with  the  mechanical  support  springs. 

By  varying  the  relative  magnitudes  of  the  persis¬ 
tent  currents  the  stiffness  of  the  two  proof  mass 
suspensions  can  be  adjusted  to  be  identical.  In 
this  case  if  the  two  experience  a  cornnon  accelera¬ 
tion  along  their  sensitive  axes  they  move  the  same 
distance  and  there  ls  no  change  in  the  inductance 
of  the  differential  pickup  coil.  A  fourth  coll  on 


SQUID 


Figure  3.  Common  Mode  Detection  Circuit 
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the  toroidal  transformer  la  provided  to  vary  the 
difference  in  the  two  peralatent  currenta .  A 
current  in  this  coil  couplea  flux  in  oppoaite 
directions  into  the  two  pickup  loops.  The  currents 
vary  to  maintain  the  net  flux  constant . 

IV.  TEST  RESULTS 

Initial  tests  have  been  performed  to  check  the 
basic  properties  of  the  instrument.  The  frequency 
of  the  differential  mode  is  62.7  Ht,  and  it  has  a 
Q  od  400.  This  Q  is  limited  by  the  construction 
of  the  gradlometer  rather  than  the  intrinsic  Q  of 
the  material.  The  Q  of  400  is  sufficient  to  reduce 
thermal  noise  below  1E/*^uT 

The  scale  factor  of  the  Instrument  was  measured 
using  the  Internal  tilting  system  with  2.0A  stored 
in  the  differential  pickup  loop.  An  experimental 
valve  of  0.30  pV/E  was  measured.  This  compares 
well  with  the  theoretical  value  of  0.39  pV/E.  The 
reduction  is  probably  due  to  flux  leakage  in  the 
coupling  transformer. 

Conanon  mode  rejection  was  studied  by  tilting 
the  cryostat  by  10  radians  and  observing  the 
response  of  the  differential  readout.  It  was  found 
that  the  response  could  be  driven  through  zero  by 
adjustment  of  the  currents  in  the  common  mode  pickup 
circuit  as  expected.  The  best  matching  obtained 
corresponded  to  a  mismatch  in  the  spring  constants 
of  5  x  10 

The  noise  power  spectrum  of  the  instrument  was 
examined  under  several  circumstances.  During  most 
experiments  the  cryostat  was  hung  on  latex  tubing 
to  form  a  low  pass  mechanical  filter  with  a  reso¬ 
nant  frequency  of  1.3  Hz.  With  this  filter  the 
noise  level  between  2  Hz  and  10  Hz  was  below 
20  E//Hz.  Below  1  Hz  the  noise  level  rises  to 
550  E/SHz.  The  sources  of  noise  have  not  been  fully 
identified;  however,  we  expect  that  the  large 
increase  in  noise  at  low  frequencies  Is  due  to 
thermal  fluctuations.  We  expect  this  to  improve 
with  better  temperature  control.  We  also  expect  to 
Improve  the  instrument's  rejection  of  seismic  noise 
by  cooling  the  system  In  low  magnetic  field.  This 
will  eliminate  pick-up  from  signal  leads  moving  with 
respect  to  trapped  flux  in  the  instrument. 


V.  INVERSE  SQUARE  UW  TEST 
Recently  it  has  been  suggested  thet  the  Inverse 
square  lew  of  gravity  may  be  violated  at  laboratory 
dlatancea  (Long,  1976).  We  are  planning  to  use  our 
gradlometer  in  a  nearly  null  experiment  to  search 
for  such  a  deviation.  Since  the  gradlometer  only 
measures  gradients  relative  to  an  arbitrary  zero, 
we  must  devise  a  ti  :  mass  geometry  which  produces 
a  nearly  zero  gravity  gradient  and  compare  the 
output  of  the  in  trument  with  and  without  the  test 
mass  present. 

The  test  mass  geometry  we  have  chosen  is  a 
cylindrical  shell  large  enough  to  fit  over  our 
cryostat.  The  gravity  gradient  at  the  center  of 
this  shell  in  the  plane  perpendicular  to  the  cylin¬ 
der  axis  is  given  by 

.  CM  .  12R./.Lj_? 

*  2  ?  3/2 

R  L  ll+(2R/L)Z)  ’ 

where  M,  R,  and  L  are  the  cylinder  mass,  radius 
and  length.  If  L>>R  this  gradient  is  nearly 
zero.  Our  cylinder  will  have  a  mass  of  3.2  x  103 
kg,  a  length  of  5.)  m  and  a  radius  of  25  cm  so  that 
L  “  20  R.  The  central  gradient  is  6.7  E.  We  hope 
to  measure  this  to  an  accuracy  of  1Z.  This  would 

4 

allow  us  to  see  a  deviation  of  1  part  in  10  from 
the  inverse  square  law. 
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A  SUPERCONDUCTING  CMV1TY  GRADIOH1TERS  -  t.  Mapoles 
DISCUSSION 

BROLLEY — Based  on  your  tup*-  .-sent,  what  1  in  t  can 
you  place  on  the  reat  m«  *  s  .»i  the  graviton’ 

HAPOLES--U  depends  on  (he  range  in  whici.  ve  do  the 
experiment.  We  can't  really  place  leith  -  an  upptr 
or  lower)  a  limit  on  che  reet  maaa.  With  this  type 
of  experiment,  we  can  find  the  ran,'*  of  the 
violation  to  be  outside  of  the  range  between  10  rm 
and  1  m.  this  corresponds  to  a  rang.'  of  reat 
masses  for  some  particle  other  than  th<  graviton. 
The  graviton  is  well  known  Co  be  assslet  ■ 

PAIR — There  are  some  indications  that  another  type 
of  graviton  exists. 

HAPOLES--The  only  reasonable  way  to  interpret  a 
violation  is  to  say  chat  there  is  another  particle 
which  has  a  reat  mass. 

PAIR — There  is  a  theory  which  aaya  there  are  scalar 
gravitons  with  mass. 

BROLLEY --Can  you  separate  spin  differences? 


MAFOLES  ns  can  inlet  something  about  the  spin 
•  coin  the  ,ign  of  the  force.  By  measuring  the  sign 
of  'lie  violation,  one  car.  discover  something  about 
the  spin. 

nVERTON--When  you  cool  'he  evsteB  down  initially, 
unless  you  cool  it  down  it:  a  Helmholtt  coil  or  in  a 
magnetic  shield,  you  will  have  froten-in  flux 
within  the  superconducting  system.  Does  this 
affect  the  pickup  from  your  coil?  Does  the  trapped 
flux  move? 

MAPOLES-'Yes.  That  is  the  leading  problem  in  the 
pickup.  In  fact,  it  is  even  worse  because  1  have 
Actuators  which  are  magnets.  They  have 
Superconducting  shields  of  their  own,  but  a  little 
more  than  the  earth’s  field  gets  trapped  in  the 
Superconducting  structure.  As  the  superconducting 
Structure  rotates  with  respect  to  those  magnets,  we 
do  get  a  spurious  effect.  1  think  we  can  easily 
eliminate  that  problem. 
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ABSTRACT 


The  inverse  square  law  of  gravitation  is  known  to  agree  with  astro¬ 
nomical  data  to  very  high  accuracy,  but  recent  theoretical  and  experi¬ 
mental  work  indicate  that  the  inverse  square  law  may  be  violated  at 
distances  less  than  1000  km.  Such  a  violation  would  signal  the  exist¬ 
ence  of  a  new  force. 

In  order  to  check  the  inverse  square  law  we  are  preparing  to  search 
for  a  non-newtonian  force  in  a  cylindrical  shell.  The  cylindrical  shell 
has  the  property  that  the  newtonian  effects  nearly  cancel,  so  that  we  are 
doing  a  nearly  null  experiment. 

We  have  developed  a  superconducting  gradiometer  to  measure  the  gravi 
tational  force  gradients  at  the  center  of  this  cylindrical  shell.  By 
measuring  both  the  vertical  and  horizontal  gradients  we  can  eliminate 
effects  due  to  imperfections  in  the  cylinder.  This  thesis  describes  the 
instrument  and  its  performance  as  well  as  calculations  of  the  sensitivity 
of  the  inverse  square  law  test. 
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INTRODUCTION 


Newton's  inverse  square  law  of  gravitation  is  known  to  agree  with 

astronomical  data  to  very  high  accuracy  and  has  been  widely  accepted  as 

the  correct  force  law  in  the  static  limit  at  all  mass  separations.  However, 

it  has  recently  been  pointed  out  that  existing  experimental  data  allows 

3 

large  violations  of  the  inverse  square  law  at  distances  less  than  10  km 
[Mikkelsou  1977],  [Long  1974],  and  several  authors  have  proposed  that  the 
gravitational  potential  <j>(r)  at  a  distance  r  from  a  point  source  of  mass 
M  has  the  form 


,  .  ,  GMr  .  .  “Sr, 

4>(r)  =  “  — l  1  +“e  ]  ,  [  I>;L] 

where  a  is  a  constant  greater  than  -1  and  8  and  G  are  constants  [Wagoner  1970], 
[Fujii  1971],  [O'Hanlon  1972].  This  potential  has  the  property  that  if  » 1 
the  effect  of  the  additional  term  is  negligible  and  the  inverse  square  law 
is  recovered  even  if  a  is  of  the  order  of  unity.  The  gravitational  force  on 
a  test  mass  m  in  this  potential  is  given  by 

f  -  G(r)>fa  [1.2] 

2 

r 

where 

G(r)  =  G  [1  +a(l  +  gr)e]  _6r  .  (  1.3] 

Thus,  if  one  measures  the  force  between  two  particles  as  a  function  of  sepa¬ 
ration  in  the  presence  of  the  potential  [i.l]  one  would  find  that  the  gravi¬ 
tational  constant  G,  as  it  is  usually  defined  assuming  an  inverse  square  force 
law,  is  actually  a  function  of  mass  separation.  Long  [Long  1976]  claims  to 
have  observed  such  a  variation  in  the  gravitational  constant.  The  expression 


for  G(r)  in  [l.3|  has  the  interesting  feature  that  for  Br  »  1,  G(r)  =  G, 
but  if  Br  <<  1,  G(r)=G(l+a).  Values  of  the  gravitational  constant  are 
obtained  from  laboratory  measurements  with  mass  separations  much  less  than 
100  m,  so  if  100  m  <  S  ^  < 10^  m;  then  the  value  of  G  to  be  used  in  astro¬ 
nomical  calculations  is  not  the  value  measured  in  the  laboratory.  Evidence 
for  and  against  a  violation  is  reviewed  in  Chapter  3. 

The  possibility  of  a  violation  of  the  inverse  square  law,  and  the  need 
for  more  accurate  laboratory  tests  has  motivated  the  work  of  this  thesis. 
Using  techniques  of  low  temperature  physics,  we  have  developed  an  instrument 
for  measuring  gravitational  gradients.  The  design  and  construction  of  this 
instrument  are  described  in  Chapter  1.  The  results  of  initial  testing  are 
described  in  Chapter  2. 

This  instrument  is  intended  for  use  in  an  experimental  test  of  the 

inverse  square  law.  This  test  consists  of  measuring  the  three  components 

_  2  ? 

of  ;  i  at  the  center  of  a  cylindrical  shell  and  verifying  that  V  <p  is  zero. 

Room  temperature  gravity  gradiometers  have  been  developed,  [Forward  1973], 

[  Bell  1977)  ,  [  Trageser  1975]  ,  but  these  instruments  are  not  able  to  check 
that  7^5  is  zero  because  they  measure  sheer  gradients  or  linear  combinations 
of  sheer  and  in-line  gradients.  The  gravitational  gradients  produced  by  an 
inverse  square  force  law  as  well  as  those  produced  by  the  Yukawa  potential 
of  equation  [  1 . 1 ]  in  a  cylindrical  shell  are  calculated  in  Chapter  4. 

Various  error  sources  are  also  considered  in  Chapter  4.  These  include  posi¬ 
tioning  errors,  errors  due  to  the  finite  size  of  the  instrument,  and  errors 

2 

due  to  imperfections  of  the  cylinder.  By  verifying  that  V  $  is  zero  one 
can  eliminate  all  errors  due  to  imperfections  of  the  cylinder.  In  fact 
one  could  use  any  source  geometry  to  check  the  inverse  square  law  if  one 
checks  that  7“5  =  0  (  Paik  1979].  However,  the  cylinder  geometry  offers 
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several  large  advantages  because  the  gradients  at  the  center  of  the 
cylinder  produced  by  an  inverse  square  force  law  are  attenuated  by  the 
ratio  of  the  diameter  to  the  length  squared.  For  our  cylinder  this 
factor  is  0.01.  Because  of  this  the  accuracy  with  which  the  instrument 
must  be  calibrated  is  also  reduced  by  100.  Another  advantage  of  the 
cylinder  is  the  uniformity  of  the  gradients  inside.  Because  the  gradients 
are  so  uniform  near  the  center,  the  accuracy  with  which  the  gradiometer 
needs  to  be  positioned  is  very  low,  and  the  effects  of  the  finite  size  of 
the  gradiometer  are  also  much  smaller  than  they  are  with  other  source 
geometries . 

The  size  of  the  gradients  in  the  cylinder,  and  the  desire  to  measure 
vertical  as  well  as  horizontal  gradients  have  influenced  the  design  of  the 
gradiometer.  The  gradiometer  consists  of  two  accelerometers  whose  sensitive 
axes  lie  along  the  same  line.  Relative  motion  of  these  masses  is  sensed 
so  that  a  differential  force  over  the  separation  of  the  masses  is  measured. 
This  differential  force  or  acceleration  is  an  approximation  of  the  true 
gradient.  The  proof  masses  are  mechanically  suspended  with  pairs  of  springs 
which  provide  very  low  compliance  in  all  directions  except  along  the  sensi¬ 
tive  axis  of  the  gradiometer.  However,  there  is  a  limit  to  how  large  a 
compliance  can  be  allowed  along  the  sensitive  axis  if  the  gradiometer  is  to 
be  operated  vertically  and  horizontally,  because  the  proof  masses  will  sag 
when  the  gradiometer  is  turned  from  horizontal  to  vertical.  The  need  to 
keep  the  support  springs  relatively  stiff  puts  a  burden  on  the  readout, 
since  the  displacements  produced  by  gradients  are  reduced  by  stiff  springs. 
For  this  reason  we  make  use  of  an  extremely  sensitive  superconducting 
readout  described  in  section  1.2. 
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To  get  some  idea  of  the  displacement  sensitivity  required,  let  us 

calculate  the  displacement  produced  by  an  acceleration  gradient  of  1  E, 

-9  2 

where  1  E  is  a  change  in  acceleration  of  10  m/s  over  a  distance  of  1  m, 
-9  2 

or  10  /s  .  The  resonant  frequencies  of  the  spring-mass  systems  which 
make  up  the  accelerometers  are  60  Hz,  and  their  separation  is  3.2  cm,  so 
a  gradient  of  IE  produces  a  differential  motion 


Ax 


(3.2  cm)(10~9s  2) 

2  2-2 
(2tt)  (60)  s 


2.3  x  10 


-14 


cm . 


This  is  an  extremely  small  displacement,  but  it  is  detectable  with  the 
inductance  modulated  readout  and  the  SQUID  magnetometer. 

With  a  60  Hz  resonant  frequency,  the  sag  of  an  accelerometer  under 
one  g  is 


Ax 


981  cm/s2 
(2tt)2(60)2/s2 


6.9  x  10 


cm . 


This  number  is  small  enough  to  allow  the  gradiometer  to  be  oriented  verti¬ 
cally  or  horizontally  without  changing  internal  components  of  the  instrument 


CHAPTER  1 


The  Superconducting  Gravity  Gradiometer 


/,  ,  %  JS 


1  *  1  The  Gravity  Gradient  Tensor 

According  to  the  classical  theory  of  gravitation,  the  gravitational 
field  in  vacuum  is  described  by  a  scalar  potential  tj>(r)  which  satisfies 
Poisson's  equation.  The  force  on  a  test  particle  of  mass  m  is 


F  =  -m  $<P  (r  )  . 


[1.1] 


The  tidal  forces  generated  by  <}>  are  described  by  the  gravity  gradient  tensor 
F,  given  by 
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The  elements 

of  r 

are  subject 

to  four 

constraints  since 
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I 
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[1.5] 

• 

• 

/  \V 

since  <J>  satisfies  Poisson  s  equation.  This  leaves  five  independent  compo¬ 


nents.  In  addition  to  gradients  due  to  the  gravitational  field,  there  are 
other  gradients  present  if  the  sensor  is  not  operated  in  an  inertial  refer¬ 
ence  frame.  These  affects  basically  arise  from  the  centripetal  forces  on 
the  proof  masses  when  the  sensor  rotates  with  respect  to  inertial  space. 
Suppose  that  the  gradiometer  is  fixed  in  some  noninertial  frame  0’  as  shown 
in  figure  1.1.  If  r  is  the  gradiometer  position  in  the  inertial  frame  0 
and  r'  is  its  position  in  the  O'  frame,  and  the  O'  frame  has  an  angular  velo 
city  co,  then 

••  •  • 

d2r/d t2  =  r'  +  2co  x  r '  +  co  x  r '  +  co  x  (co  x  r'). 

[1.6] 

If  the  gradiometer  is  at  rest  in  O'  the  first  two  terms  are  zero,  then 


-►2  .  2  -► 

dr/dt  =  ai  x  r'  +  0)  x  (a)  x  r  1 )  . 


[1.7] 


To  calculate  the  acceleration  gradients  we  need  to  calculate 

,2-*- 


ij  3r'.  [_dt2  •  J  J 


The  vector  co  is  fixed  in  the  primed  frame  so 


cox  r  '  =  i  '  (co  z '  -  coy')  -  j  '  (co  z  '  -  co  x  ' ) 
y  z  x  z 


+  k(co  y'  -  cox'), 
x  y 


and 


to  x  (co  x  r ')  =  i '  (co  +co  )x'  +  co  oj  y'  +  u  co  z'l 

L  y  z  y  x  z  j 


co  co  z 
y  z 


'] 


+  jf  rajcoxf-(a)  -f  co  ;yf  + 
y  x  x  z 

+  k'  fco  co  x’  +  co  co  y'  -  {u  2  +  co  2)z'l  . 


1.8] 


[1.9] 


[  1.10] 


Taking  the  partials  we  have 
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[1.11] 


The  instrument  which  we  will  describe  here  measures  a  diagonal  component 

2  2 

of  T  or  I”,  hence  we  are  interested  in  error  terms  of  the  form  co^  +  co^  . 
Suppose  that  the  sensitive  axis  is  along  the  x  direction  and  the  instrument 
is  given  a  tilt  about  the  z  axis 


Then 


0=0  coscot  . 
o 


4)  =  8  =  0)0  sin0)t, 

z  o 


[1.12] 


[1.13] 


and 


e  V 


(1  -  cos2o)t)  , 


[1.14] 


2  2 

The  gradient  consists  of  a  component  at  2o)  and  a  d.c.  component  0q  o)  /2. 

This  later  term  represents  a  serious  source  of  error  since  it  increases  as 

the  square  of  the  frequency.  For  instance,  an  angular  velocity  of  3  xlO 

-9 

rad/sec  or  6°  per  hour  produces  an  equivalent  gravity  gradient  IE  =  10  /sec 


1 . 2  The  Superconducting  Readout 

The  superconducting  readout  circuit  is  shown  in  figure  1.2.  It  consists 
of  an  inductor  rigidly  attached  to  one  of  the  proof  masses  and  coupled  to 
a  SQUID  by  a  superconducting  transformer.  is  located  as  closely  as 
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SQUID 


Figure  1.2  Differential  Readout  Circuit 


possible  to  the  surface  of  the  other  proof  mass  so  that  its  inductance  is 
strongly  modulated  by  relative  motion  of  the  two  proof  masses.  A  persis¬ 
tent  current  is  stored  in  the  loop  consisting  of  and  L^.  Since  the 
total  flux  in  this  loop  must  remain  constant,  a  change  in  the  inductance 
of  causes  a  change  in  the  current  which  is  detected  by  the  SQUID. 

The  inductor  is  a  spiral  of  niobium  wire  wound  on  a  flat  dielectric 
substrate.  This  flat  circular  coil  is  glued  to  one  of  the  proof  masses. 

A  flat  superconducting  plane  which  is  attached  to  the  second  proof  mass  is 
located  about  5  mills  from  the  coil.  The  readout  coil  is  shown  in  figure 
1.7.  According  to  Ampere's  law  the  field  B  produced  by  a  current  I  in  the 
coil  is 

B  =  n|i  I,  [1.15] 

o 

where  n  is  the  number  of  turns  per  unit  length  in  the  coil.  This  field 
occupies  a  volume  As  where  A  is  the  area  of  the  coil  and  s  is  coil  proof 
mass  spacing,  so  that 

L,  =  p  n2  As.  [  1.16] 

1  o 


This  can  be  rewritten  as 


where 


2AJ 
M  n  Ad 
o 


[1.17] 


[1.18] 


and  d  is  the  equilibrium  spacing  between  the  proof  masses. 


We  can  write  the  equations  of  motion  of  the  two  proof  masses  as 


mx^  +  kx^  =  -  a  I  /2 


mx2  +  kx2  =  a  I  /2, 


fl.19] 


where  a  =  L^°/d  =  U^n^A.  We  have  ignored  any  damping  and  any  mechanical 
mismatch  between  the  proof  masses.  If  £  =  x^  -  x^  we  have 


»*“_  •  w*  » 


v  .  2  al“ 

C  +  co  5  =  — 
o  m 


[1.20] 


where  a)  =  k/m.  As  current  is  stored  in  L,  the  magnetic  field  pushes  the 
o  1 

proof  masses  apart  and  d  increases.  If  h  is  the  spacing  with  zero  current. 


,  ,al 

h  +  o 


[1.21] 


where  Iq  is  the  mean  value  of  the  persistent  current  and  d  is  the  mean  value 
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of  the  coil  proof  mass  spacing.  Now  if  we  define  £  =  £  -  al^  /njW0  the 


equation  of  motion  becomes 


5+  ao‘C 


[l2  -  ’o2]  • 


[1.22] 


Since  flux  in  any  superconducting  loop  is  fixed,  we  have  two  constraints 


*1  =  (Lx  +  L2)I  -  Mi  =  (Ll  +  L2)Iq 


[ 1.23a] 


^  =  (L,  +  L.)i  -  MI  =  -MI 

2  3  4  o 


[1.23b] 


where  M  is  the  mutual  inductance  of  the  transformer  in  figure  1.2.  The  mean 
value  of  i  is  forced  to  be  zero  when  current  is  stored  with  a  heat  switch 


which  causes  part  of  the  loop  to  go  normal.  Solving  for  I  and  i,  we  find 
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.  .< 


io  1  _ 

1(4  +  L2)(L3  +  4)  -  M2  ] 


"V^l  -  L1  > 

(4  +  L2)(L3  +  4)  - 


[1.24] 


[1.25] 


If  we  let  (4°  +  4X4  +  L4)  -M2  =  L2  and  expand  4  =  4°|l  +  j-J,  then 


to  lowest  order  in  5 


'“'I  L  « 


[1.26] 


Using  [1.24]  we  can  write  the  equation  of  motion  [1.22]  as 
",  2  „  a  2  [  a(L3+L4)  c]  oil  2 

5  +  Wo  «  =  T  ro  l1  +  - J—  »  "  -2-  • 


[1.27] 


Expanding  the  right  side  and  keeping  the  lowest  order  in  C,  we  have 


,  -2a2I  2(L_+L  )£ 

r  t  2  r  o  3  4 

S  +  CO  5  - - o - 

0  mL2 


[1.28] 


This  expression  gives  the  frequency  shift  of  the  differential  mode  with 
persistent  current  I  ,  but  L  is  not  independent  of  I  .  So  let  us  write 


and  define 


T  0  ,  ,  o 

4  =  a  h  +  - — 

mo) 

o 


2  2 
L  =  (ah  +  L_)(L-+L.)  -  M 

o  Z  O  4 


[1.29] 


[1.30] 


Then  we  can  write  [1.28]  as 


£  +  CD 


-2“V(L3  +  L4>  (,  ,  °V<W  ) 

i.  =  - X -  U  +  - o - X -  I 
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or  we  can  write  frequency  of  the  differential  mode  f  as 
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2  2  2 
The  behavior  of  f  versus  I  is  shown  in  figure  1.3.  f  increases 

2  2  2  2 

rapidly  initially  with  a  slope  of  2  a  I  (L.  +  L,)/(2tt)  mL  but  quickly 

o  J  4  o  ‘ 

2 

rolls  off  and  approaches  the  asymptotic  value  of  3f  .  To  find  an  optimal 

o 

value  of  I  ,  note  that  a  gravity  gradient  G  produces  a  differential  displace- 

2 

ment  of  the  proof  masses  Gb/(2rrf)  where  b  is  the  gradiometer  baseline.  So 


from  [  1 . 26] 
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From  this  it  follows  that  i/G  is  maximized  when 
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Substituting  this  into  the  brackets  of  [1.33]  gives 
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where  Lc  =  ah.  To  impedance  match,  we  construct  the  transformer  so  that 
L2/Lc  =  L^/L^  =  y,  and  we  assume  that  M  ~  /L^L^.  Then 


( 1+  2y )  h  2oj 


i 2  nr 

o _  I  __c_ 

h  2.  2  V  L4 


[1.36] 


is  the  scale  factor  of  the  gradiometer.  This  is  optimized  for  Y>>1» 

but  any  value  greater  than  one  gives  a  scale  factor  close  to  the  optimal. 

One  consequence  of  [1.36]  is  that  the  scale  factor  increases  as  0)q  ^ 

_2 

rather  than  oo  if  the  resonant  frequency  of  the  gradiometer  proof  masses 
o 

is  reduced  keeping  their  masses  constant,  since  I  is  proportional  to  . 
This  is  a  consequence  of  the  current  causing  the  gap  between  the  proof  masses 
to  increase  as  well  as  the  increase  of  spring  constant  provided  by  the  elec¬ 
trical  spring.  Another  interesting  consequence  of  [ 1 . 36]  is  that  the  scale 
factor  increases  only  as  1//  h  since  is  proportional  to  h. 

1 . 3  Description  of  the  New  Gradiometer 

Our  primary  goal  in  the  design  and  fabrication  of  our  gravity  gradio¬ 
meter  was  to  minimize  the  instrument's  sensitivity  to  linear  accelerations 
both  along  the  proof  mass  direction  and  perpendicular  to  it.  Linear  accel¬ 
erations  in  the  laboratory  exert  forces  on  the  proof  masses  which  are  very 
large  compared  to  the  gravitational  forces  we  would  like  to  measure,  so 
careful  balancing  of  the  gradiometer  is  needed  to  eliminate  them  as  a 
noise  source. 

In  addition,  we  decided  to  build  an  in-line  gradiometer  which  measures 

one  of  the  diagonal  components  of  T.  As  we  will  show  in  Chapter  A,  it  is 

2 

highly  desirable  to  be  able  to  measure  all  three  components  of  V  <j>  so  we 
chose  a  method  for  suspending  the  proof  masses  which  allows  us  to  operate 
the  gradiometer  with  its  sensitive  axis  vertical  as  well  as  horizontal. 

In  order  to  make  the  balancing  problem  manageable  we  chose  a  mechanical 
design  which  enabled  us  to  use  only  one  coil  in  the  readout  system  described 
in  section  1.2.  Experiments  have  been  done  with  readouts  using  several 
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coils,  [  Paik  1978]  but  in  these  gradiometers  careful  matching  of  the  coils 
as  well  as  the  proof  masses  is  necessary  to  reject  linear  accelerations. 

Basically  the  gradiometer  consists  of  two  superconducting  proof  masses 
adjacent  to  each  other.  A  single  niobium  coil  is  mounted  in  the  gap  between 
the  proof  masses  so  that  the  coil  inductance  is  proportional  to  their  spacing. 


»  « 


This  readout  system  is  analyzed  in  detail  in  the  previous  section.  Since 
the  sensitive  axes  of  the  proof  masses  are  parallel,  the  gradiometer  detects 
changes  in  one  of  the  diagonal  components  of  the  gravity  gradient  tensor. 


Each  of  the  two  proof  masses  is  supported  with  two  mechanical  springs. 
The  mechanical  springs  are  folded  cantilevers  cut  into  circular  disks  of 
niobium.  A  folded  cantilever  with  two-fold  symmetry  is  shown  in  figure  1.4. 
In  this  case  the  thickness  of  the  plate  is  reduced  in  the  two  rectangular 
areas  hatched  in  the  figure  in  order  to  restrict  the  bending  to  those  areas. 
Slots  are  Chen  cut  in  the  indicated  locations.  These  slots  allow  the  center 
to  move  perpendicular  to  the  plane  of  the  plate  as  shown  in  the  lower  part 
of  the  figure.  Since  the  center  of  the  spring  is  free  to  move  in  the  plane 
of  the  plate,  the  total  length  of  the  spring  remains  constant  and  the 
restoring  force  is  provided  by  pure  bending  without  stretching  of  the  spring 
material.  The  stiffness  can  be  calculated  using  the  cantilever  equation 
[  Landau  1959]  .  In  the  case  shown  the  spring  displacement  5  is 


wEh  [1.37] 

where  f  is  the  force  in  the  center,  h  is  the  cantilever  thickness,  E  is 
the  Young's  modulus  of  the  material  and  Z  and  w  are  the  lengths  shown  in 
figure  1.5.  The  spring  constant  k  is 


m 
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The  springs  used  in  the  gradiometer  have  three-fold  symmetry  rather 

than  two  and  a  picture  of  one  of  them  is  shown  in  figure  1.5.  Since  two 

such  springs  are  used  to  support  each  proof  mass,  there  are  three  pairs  of 

cantilevers,  and  the  net  spring  constant  for  each  proof  mass  is 

,  3  Ewh3 

k  =  7  — - —  ■ 

4  S.  [1.39] 

These  springs  were  fabricated  from  solid  disks  of  niobium  10.5  cm  in 
diameter  and  0.32  cm  thick.  The  faces  were  carefully  machined  to  produce 
a  final  thickness  uniformity  of  approximately  50  microns.  The  rectangular 
depressions  were  milled  to  within  a  few  thousandths  of  an  inch  of  their 
final  thickness,  and  the  part  was  then  stress  relieved.  In  the  next  machining 
stage,  the  rectangular  depressions  were  cut  to  their  final  thickness  of 
0.071  cm  using  electron  discharge  machining  (EDM)  with  a  rectangular  elec¬ 
trode  of  copper  and  graphite.  EDM  was  used  to  achieve  a  careful  match  of 
all  of  the  cantilevers,  and  thereby  match  the  final  proof  mass  frequencies. 

In  the  final  stage  of  machining  the  narrow  slots  were  cut  using  wire  EDM. 

In  this  process  wire  is  continuously  circulated  through  the  part  and  used 
as  an  electrode.  The  part  is  moved  on  a  numerically  controlled  table  as 
the  slot  is  cut. 

The  proof  mass  subassembly  consists  of  two  springs,  two  cover  pieces, 
an  annular  mass,  and  a  ring  used  to  hold  the  edges  of  the  springs  apart. 

The  cover  pieces  are  circular  disks  of  niobium  8.51  cm  in  diameter  which 
are  flat  on  one  side  and  have  a  threaded  stud  on  the  opposite  side.  The 
flat  sides  are  designed  to  modulate  the  inductance  of  the  readout  coil  or 
tuning  coils  described  in  section  1.2.  The  flats  were  polished  to  ensure 


flatness  and  isotropy  so  that  close  spacing  to  the  coils  could  be  obtained, 
and  so  that  the  inductance  of  the  coils  is  modulated  only  by  linear  motion 
of  the  surfaces.  The  threaded  studs  screw  into  an  annular  disk  of  niobium 
which  provides  the  bulk  of  the  mass  for  each  proof  mass.  A  spring  is  captured 
between  each  cover  piece  and  proof  mass.  An  exploded  view  of  one  of  the  two 
proof  mass  subassemblies  is  shown  in  figure  1.6.  The  ring  which  is  used  to 
maintain  the  spacing  of  the  edges  of  the  springs  is  thicker  than  the  proof 
mass  so  that  each  spring  is  biased  away  from  its  equilibrium  position  and 
never  passes  through  its  zero  point.  This  eliminates  certain  forms  of  non¬ 
linear  behavior.  Since  the  springs  have  opposite  biases,  the  first  order 
nonlinear  terms  also  cancel.  The  ring  was  ground  to  its  final  thickness 
to  keep  its  two  sides  parallel  to  one  part  in  10^. 

By  using  a  design  in  which  all  of  the  parts  stack  up  on  each  other, 
and  all  parts  have  highly  parallel  surfaces,  we  ensure  that  the  sensitive 
axes  of  the  two  proof  masses  are  also  parallel.  This  parallelism  is  main¬ 
tained  in  all  the  parts  to  a  few  parts  in  10^ . 

The  readout  coil  is  shown  in  figure  1.7.  It  is  a  single  layer  of 
0.0089  cm  diameter  niobium  wire  wound  on  a  0.25  cm  thick  disk  of  macor. 

The  coil  has  400  turns  and  a  diameter  of  6.9  cm.  The  lead  to  the  center  of 
the  coil  is  glued  into  a  slot  in  the  macor  substrate  so  that  both  leads 
come  off  the  edge  of  the  coils.  The  macor  coil  form  is  glued  to  one  of 
the  four  cover  pieces.  This  cover  piece  is  0.06  cm  thinner  than  the  other 
three  to  compensate  for  the  coil  form  mass. 

Once  the  two  proof  mass  subassemblies  were  constructed,  they  were 
stacked  on  top  of  each  other  with  a  ring  holding  them  apart.  The  thickness 
of  the  ring  was  chosen  so  that  the  readout  coil  is  about  0.013  cm  from  the 
opposite  proof  mass.  This  structure  was  mounted  on  an  aluminum  fixture  which 
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Figure  1.9  Gradiometer  Cross  Section 
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maintained  the  angular  orientation  of  the  parts.  A  niobium  housing  tube 

9.3  cm  long  with  a  10.5  cm  inside  diameter  was  placed  over  the  parts. 

The  ends  of  the  assembly  are  sealed  with  a  pair  of  tuning  and  feedback 
coils  shown  in  figure  1.8.  The  inner  coil  on  each  form  is  currently  used 
to  match  the  resonant  frequencies  of  the  proof  mass  subassemblies  as 
described  in  section  1.5.  The  annular  outer  coil  was  included  for  future 
use  as  a  feedback  coil.  Fixtures  for  making  superconducting  joints  and 
heat  switches  are  mounted  on  the  opposite  sides  of  these  tuning  coils. 

There  is  also  a  superconducting  transformer  mounted  on  each  tuning  coil; 
one  for  the  differential  readout  circuit  and  one  for  the  common  mode  readout 
circuit . 

The  entire  assembly  is  locked  together  with  a  pair  of  rings  which  screw 
into  either  side  of  the  housing  tube.  A  cross  section  of  the  assembled 
gradiometer  is  shown  in  figure  1.9.  The  assembled  gradiometer  has  a  baseline 
of  3.2  cm. 

All  of  the  gradiometer  components,  except  for  the  dielectric  coil  forms, 

-3 

are  niobium.  Since  niobium  undergoes  a  contraction  of  1.4  x  10  when  cooled 
to  4.2°K,  all  parts  must  be  made  of  niobium  to  prevent  differential  contrac¬ 
tion  from  altering  the  geometry. 

1 . 4  The  Gradiometer  Support  Assembly 

The  gradiometer  is  supported  on  an  aluminum  and  titanium  fixture  in  the 
cryostat  which  allows  the  gradiometer  axis  to  be  positioned  vertically  or 
horizontally.  This  fixture  holding  the  housing  tube  is  shown  in  figure 
1.10.  The  housing  tube  is  clamped  between  two  titanium  pieces  which  bolt 
together  above  and  below  the  gradiometer  as  in  figure  1.11b.  This  part  is 


clamped  to  the  bottom  of  a  pair  of  cross  flex  hinges  shown  in  1.11a.  These 
hinges  allow  the  gradiometer  to  be  tilted  inside  the  vacuum  chamber  of  the 
cryostat  to  produce  test  signals.  Titanium  is  used  for  the  parts  which 
make  contact  with  the  gradiometer  because  its  thermal  contraction  is  a 
good  match  to  that  of  niobium. 

Mechanical  force  is  transmitted  to  these  hinges  by  a  pair  of  arms 
inside  the  aluminum  support  beams.  These  arms:  are  rigidly  attached  to  the 
gradiometer  with  an  arm  that  a  mechanical  force  is  applied  to  in  order  to 
apply  torque  to  the  hinges. 

Two  methods  of  applying  torque  to  the  assembly  have  been  used.  In 
the  f  irst  experiments  with  the  gradiometer  a  pair  of  superconducting  sole¬ 
noids  were  attached  to  the  torquer  arm  and  inserted  into  the  gaps  of  a  pair 
of  permanent  magnets  fixed  to  the  top  plate.  Current  with  opposite  polarity 
in  the  solenoids  produces  an  angular  deflection  of  the  gradiometer.  This 
arrangement  was  replaced  with  a  pair  of  PZT  bimorphs  later  because  flux 
from  the  magnets  was  trapped  in  the  gradiometer  as  it  was  cooled  through 
its  transition  temperature,  and  this  flux  may  have  caused  additional 
coupling  of  the  readout  to  mechanical  motion. 

1 . 5  Magnetic  Tuning  and  Common  Mode  Detection 

Even  with  excellent  machining  tolerances  the  mechanical  springs  can 
only  be  matched  to  a  few  percent.  Since  we  have  designed  the  system  to 
reject  cross  accelerations  to  a  few  parts  in  10^  we  would  like  to  improve 
the  frequency  matching  to  this  level  also.  In  this  case  an  equal  force  on 
both  masses  will  produce  the  same  displacement,  and  the  inductance  of  the 
readout  coil  will  be  unchanged.  To  accomplish  this,  we  have  a  magnetic 


tuning  circuit  in  the  gradiometer  which  adjusts  the  relative  stiffnesses 


of  the  proof  mass  suspensions  with  respect  to  the  case. 

The  circuit  is  shown  in  figure  1.12.  The  two  inner  coils  of  figure 
1.8  are  attached  to  matched  inductors  of  a  toroidal  core.  A  third  winding 
on  this  core  is  used  to  adjust  the  currents  in  these  loops.  A  fourth 
winding  on  this  core  couples  to  a  SQUID  which  detects  common  mode  motion  of 
the  proof  masses.  Its  inductance  is  much  smaller  than  the  other  three  and 
will  be  ignored  in  calculating  the  frequency  shifts  caused  by  currents  in 
the  other  three  loops. 

The  inductances  and  are  given  by 


and 


[1.40a] 


[  1.40b] 


where  L^°  ~  L 2°  and 


df 


The  flux  equations  for  the  three  loops  are 


♦  j  =  (4+L)^  -  LI2  -  , 

<P2  =  (L2  +  L)I2  -  LI1  +  MtLt  , 

and 

<f>3  =  LtIt  +  MT(I1  -  I2)  . 


[ 1.41a] 
[1.41b] 

[1.41c] 


We  have  assumed  that  the  mutual  inductance  between  loops  one  and  two  is 


just  L,  and  we  will  also  assume  that  M^,  =  /L^L  .  Solving  for  the  currents 
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Similarly, 
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Using  the  currents  I  and  I2  we  can  write  the  equations  of  motion  for  the 
proof  masses. 
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If  we  expand  and  we  then  have 
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Hence  the  proof  mass  resonant  frequencies  are  given  by 

/.  «r  ,  \2 
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[  1.48b] 


With  current  in  the  tuning  loops,  the  electrical  springs  act  in  parallel 
with  the  mechanical  springs.  In  order  to  tune  the  gradiometer  the  sum  of 
these  two  spring  constants  must  be  equal  for  the  two  proof  masses.  There 
are  two  levels  of  adjustment  available.  By  storing  different  <j>^  and 
$2  the  two  suspensions  can  be  roughly  matched,  then  can  be  adjusted 
for  a  more  accurate  match.  When  current  is  stored  in  the  readout  loop 
the  accelerometers  are  coupled  and  the  normal  modes  of  the  system  split 
into  a  mode  which  is  primarily  differential  motion  whose  frequency  increases 
with  readout  current  as  in  figure  1.4,  and  a  mode  which  is  primarily  equal 
motion  of  the  masses.  Bui  the  tuning  problem  remains  the  same.  In  order 
to  make  the  gradiometer  Insensitive  to  common  forces  on  the  proof  masses, 
each  proof  mass  must  have  the  same  spring  constant  between  itself  and  the 
case  of  the  instrument,  so  tuning  is  accomplished  in  the  same  way. 

Detection  of  the  common  mode  is  carried  out  with  a  circuit  quite 
similar  to  the  fine  tuning  loop.  Let  us  consider  the  circuit  consisting 
of  the  two  tuning  loops  and  the  loop  which  couples  to  the  SQUID,  ignoring 
the  fine  tuning  loop  since  currents  stored  in  this  loop  are  small.  In  this 
case  the  equations  of  motion  are  very  similar  to  [  1 . 41]  , 


*1  =  (L1  +  L)I1-  LI2  'Vs  * 
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f  1.49c] 


where  Mg  =  /L^L.  The  equilibrium  value  of  Ig  is  kept  at  zero  with  a  heat 

switch.  Solving  for  I  we  find 
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Now  in  ideal  circumstances  the  balanced  gradiometer  would  have  <J>g  =  0,  = 

<f> 2 »  Lx°  ~  ^2°.  and  d^  =  d2>  Then  Ig  is  reduced  to 
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to  lowest  order  in  x^/d^  and  x^/d^.  In  the  ideal  case  Ig  responds  only  to 
common  mode  motion  of  the  proof  masses.  In  practice  <j>^  #  (j>2,  d^  ^  d^,  and 
U^°  t  1^2°  so  does  couple  to  differential  motion.  This  is  not  really  a 
problem  because  common  mode  motion  is  much  larger  than  the  differential 
motion,  and  the  signal  from  the  differential  SQUID  can  be  used  to  identify 
that  part  of  the  signal  from  the  common  mode  SQUID  corresponding  to  differ¬ 
ential  motion. 

The  common  mode  readout  circuit  was  included  in  the  instrument 
because  in  the  future  we  will  use  its  output  in  a  feedback  loop  to  stabilize 
the  positions  of  the  proof  masses.  The  feedback  signal  will  go  to  the  outer 
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coils  in  figure  1.9.  This  feedback  loop  will  reduce  coupling  of  vibra¬ 
tion  to  the  differential  readout  due  to  nonlinearity  of  the  support  springs. 
1 . 6  Feedback  Damping  of  the  Differential  Mode 

The  high  Q  of  the  differential  mode  is  a  disadvantage  for  several 
reasons.  It  causes  poor  transient  behavior  since  the  system  rings,  and 
the  large  signal  at  resonance  taxes  the  slew  rate  of  the  SQUID  detector. 

When  the  slew  rate  of  the  SQUID  is  exceeded,  it  unlocks  and  information  at 
low  frequencies  is  lost.  Consequently ,  it  is  important  to  be  able  to  damp 
this  mode.  The  most  accessible  point  at  which  to  insert  damping  is  in  the 
loop  between  the  SQUID  and  the  coupling  transformer  as  shown  in  figure  1.13. 
First  we  analyze  a  circuit  which  uses  a  real  resistance  R  in  the  SQUID  input 
loop  to  cause  loss  and  damp  the  differential  mode.  Then  a  simple  feedback 
circuit  which  synthesizes  this  resistance  is  described. 

The  equations  of  motion  for  this  system  are 
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(-L3  +  L4)  dt 


[ 1.52c] 


I  < 


where  the  symbols  are  those  used  in  section  1.1,  and  a  resistance  R  has  been 


included  as  shown  in  figure  1.13.  To  solve  for  the  frequency  we  assume 
that  £  =  £  I  =  I  +  I'e^*",  and  j  =  j't'^1'.  Then  to  first  order  in 


§  4 


£  ,  I ’  ,  and  i 1  we  have 
o 


,  2  .  2.  r  2L.°I  I’ 

( -ti)  +  OJ  )  c,  -  1  O 

o  o  - - - 

dm 


[1.53a] 
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and 


iO)  -  0)  oj^  -  +  ^3^2  =  0  ' 


[1.59] 


The  characteristic  equation  is  now  a  cubic  since  there  is  an  additional 
real  solution  which  represents  current  j  decaying  through  R  without  oscilla 

tions  in  £•  The  equation  is  easily  solved  in  two  limits.  If  0  then 

2  2  2  2 
'ji  =  ,  which  is  the  result  of  section  1.2.  If  we  have  co  =  . 

This  is  the  frequency  we  would  measure  if  the  secondary  loop  of  the  trans¬ 
former  were  removed.  In  both  cases  there  is  no  damping.  As  R  is  varied 
between  0  and  °°  the  frequency  moves  between  these  two  points,  and  the  Q 
falls  to  some  finite  value  and  returns  to  infinity. 

We  can  parameterize  [1.59]  if  we  let  x  =  ,  A  =  and  B  = 

Then 


3  2  2 

x  -  x  A  +  x  -  AB  =0 


[1.60] 


If  x^  is  a  solution,  the  Q  is  then  given  by 

I  (x  ) 

0  =  _ S _ 2 _ 

y  2Re (x  ) 
o 


[1.61] 


Equation  [1.60]  is  a  cubic  and  can  therefore  be  solved  exactly  for  the 
Q  for  any  A  and  B.  Figure  1.14  shows  the  variation  of  the  Q  with  A  for  a 
given  value  of  B  (0.98).  Figure  1.15  shows  how  the  minimum  Q  varies  with 


B.  The  result  is  approximately  that 


L  .62 


gives  the  minimum  Q. 


This  circuit  cannot  be  used  in  practice  since  the  real  solution  implies 
that  there  is  a  low  frequency  cutoff  in  the  circuit  response,  and  the  addi¬ 
tion  of  a  real  resistor  would  increase  the  noise  in  the  circuit.  However, 
the  equivalent  of  a  resistor  can  be  built  using  feedback  techniques.  This 
scheme  is  shown  in  figure  1.16.  It  consists  of  putting  the  output  of  the 
SQUID  through  a  filtering  and  phase  shifting  network  H(oj),  and  then  feeding 
back  a  current  proportional  to  the  resulting  signal  to  a  transformer  in  front 
of  the  SQUID.  This  technique  was  demonstrated  by  Kai  Wang  [  Wang  1979] ,  who 
used  this  circuit  to  reduce  the  Q  of  a  superconducting  accelerometer. 

The  transfer  function  H(w)  consists  of  a  bandpass  filter  and  phase 

shifter.  The  SQUID  control  unit  converts  the  SQUID  input  current  j  to  a 

voltage  with  a  conversion  factor  Gg.  This  voltage  is  multiplied  by  H(oi), 

and  converted  to  a  feedback  current  i^  in  the  voltage  to  current  converter 

with  a  scale  factor  G  .  So 

v 


lc  =  G  H(<u)  G  j  ’  eia)t 
f  v  s 


[1.64] 


The  voltage  across  the  feedback  transformer  is 

di 


v  =  —  M 

f  f 


f 

dt 
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So  the  effective  resistance  is  just 


Kff  =  7T  =  -MX  G  iwH(w). 


[1.65] 


[1.66] 


The  bandpass  filter  in  H((o)  has  zero  phase  shift  at  resonance  so  an 


additional  phase  shift  of  90  is  needed  to  make  real  and  positive. 

So  if 


H(a>)  = 


-C0f0) 


*  2  ’ 

+  oj 


[1.67] 


then  R  ,r  becomes 
ef  f 


M.G  G  icj,w 
f  v  s  f 
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[1.68] 


In  practice  0)^  is  set  at  the  resonant  frequency  of  the  mode  to  be  damped 


and  R  reduces  to 
eff 


R  -  MrG  G  Q.gj,. 
eff  f  v  s  f  f 


[1.69] 


at  the  feedback  frequency.  By  adjusting  G v  and  Qf ,  Rgff  can  be  adjusted  to 


reach  a  minimum  Q. 


The  first  tests  which  will  be  described  involve  the  tuning  of  the 


differential  mode,  and  the  balancing  of  the  proof  mass  spring  constants  to 
reject  linear  accelerations. 

Figure  2.1  shows  the  frequency  of  the  differential  modes  as  a  function 
of  current  between  3.0  and  6.5  amperes.  As  the  current  was  reduced  to 
values  below  3  amperes  the  differential  mode  disappeared.  This  seems  to  be 
due  to  an  interference  between  the  two  proof  masses  which  is  eliminated  as 
the  proof  masses  are  pushed  apart.  It  might  be  some  small  piece  of  material 
which  connects  the  two  masses  but  loses  contact  with  one  of  them  as  current 
is  stored.  As  current  is  stored  in  the  differential  pickup  coil,  the 
common  mode  frequency  remains  fixed  at  59.6  He. 

The  data  of  figure  2.1  can  be  fitted  to  equation  [  1.32]  to  give  values 
for  the  frequency  of  the  differential  mode  with  Iq  =  0  and  a  parameter  in 
equation  [  1 .32]  , 
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L1o2a3  +  V 

(2Tr)2mL  2d2 
o 


[2.1] 


A  fit  of  the  data  to  equation  [1.37]  gives  a  frequency 


f  =  59  Hz  , 
o 

and 

y  =  70  Hz2/A2  . 

As  described  in  section  1.3,  the  design  value  for  f  is  60  Hz  so  the  data 
are  consistent  with  the  theory  for  the  folded  cantilever  springs.  The 
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Figure  2.1  Differential  mode  frequency  shift 
as  a  function  of  persistent  current. 


observed  frequency  of  the  common  mode  of  59.6  Hz  is  also  consistent.  The 
design  values  of  the  parameters  for  y  are 


=  .79  H/m, 


.-6 


L3  +  =  9  x  10  H, 


1.07  kg, 


and 


-9  2 

1.8  x  10  H  . 


The  value  for  assumes  that  the  coupling  of  the  transformer  in  the  readout 


circuit  is  unity,  so  in  practice  we  should  expect  Lq  to  be  somewhat  larger. 


See  equation  [1.30] .  These  numbers  give  a  value 


y  =  74  Hz2/A2, 


so  again  the  experiment  is  quite  consistent  with  calculated  results. 

Looking  back  at  equation  [1.34]  ,  the  optimum  current  is  related  to  y  by 
„  ,  2 


op 


o 

3Y 


[  2.2] 


so  that  the  experimental  numbers  give 


op 


4.2  A 


This  current  corresponds  to  a  field  of  525  gauss  or  not  much  more  than  half 
of  the  critical  field  of  1000  gauss  for  niobium.  This  is  a  consequence  of 
a  conservative  choice  of  pickup  coil  parameters.  A  smaller  coil  could  have 


been  used,  thereby  increasing  I  with  little  loss  of  sensitivity. 


In  order  to  balance  the  system,  the  response  of  the  differential  SQUID 
to  a  mechanical  excitation  of  the  gradiometer  was  measured  as  a  function  of 


currents  in  the  tuning  loops  described  in  section  1.5.  Mechanical  motion 


was  generated  in  two  ways.  The  internal  PZT  actuators  described  in  section 
1.4  were  driven  at  4.0  Hz,  and  the  entire  dewar  was  driven  externally  at 
4.0  Hz  using  a  mechanical  shaker.  The  dewar  was  hung  from  the  ceiling  with 
latex  tubing  to  provide  vibration  isolation.  This  gives  a  vertical  resonant 
frequency  of  approximately  1.2  Hz  for  the  dewar  so  it  appears  as  a  free  mass 
when  driven  at  4  Hz.  Figure  2.2  shows  the  amplitude  and  phase  response  of 
the  differential  SQUID  to  a  constant  drive  signal  on  the  PZT  as  current 
in  one  tuning  loop  is  changed.  A  constant  current  of  1.0  amp  was  left  in 
the  other  loop  during  the  test. 

The  principle  feature  of  figure  2.2  is  that  at  a  specific  current  the 
amplitude  response  passes  through  zero,  and  the  phase  changes  by  180°  through 
this  point.  This  tuning  point  is  the  point  at  which  one  normal  mode  of  the 
coupled  proof  masses  is  pure  differential  motion  and  the  other  normal  mode 
is  pure  common  motion.  At  this  point  the  two  proof  masses  deflect  by  exactly 
the  same  distance  under  a  static  load.  On  either  side  of  the  tuning  point 
one  proof  mass  moves  more  than  the  other,  and  there  is  a  differential  dis¬ 
placement  for  a  linear  acceleration.  As  the  current  is  moved  across  the 
tuning  point  the  stiffer  suspension  becomes  the  weaker,  and  the  sign  of  the 
reponse  to  a  common  acceleration  changes;  hence  the  180  degree  pliase  shift. 

These  results  prove  that  the  normal  mode  structure  of  the  instrument 
is  understood,  and  that  the  balancing  control  described  in  section  1.5  is 
working.  This  is  essential  if  we  are  to  make  accurate  measurements  of 


gravitational  gradients. 


BE 


The  most  important  tests  of  the  gradiometer  are  the  measurement  of  its 
scale  factor  and  the  detection  of  an  actual  gravitational  gradient.  This 
section  describes  such  measurements  and  gives  some  noise  power  spectra  from 
which  a  reasonable  assessment  of  the  instrument's  present  capabilities  can 
be  made . 


The  first  calibration  makes  use  of  the  gradiometer 1 s  sensitivity  to 
tilts  described  by  equation  [1.14] .  If  the  gradiometer  is  oriented  along 
the  2  axis,  and  is  given  a  tilt  of  amplitude  8^  and  frequency  co  about  the 


x  or  y  axis,  a  gradient 


g(t)  = 


0  V 

o  cos  2  cot 


[  2.3] 


results.  This  is  a  very  attractive  situation  because  the  gradient  signal 
can  be  calculated  entirely  from  measured  quantities.  It  is  independent  of 
the  baseline  of  the  instrument,  and  the  location  of  the  rotation  axis. 

The  dewar  was  hung  on  latex  tubing  as  described  in  the  previous  section. 
The  mechanical  shaker  was  attached  to  the  top  of  the  dewar  with  its  driving 
axis  parallel  to  the  gradiometer  axis.  Since  the  top  of  the  dewar  is  about 
two  feet  above  the  center  of  mass  of  the  hanging  dewar,  the  a.c.  force  from 
the  shaker  applies  a  torque  to  the  dewar  and  a  small  angular  motion  results. 
This  angular  motion  is  detected  with  an  automatic  autocollimator  [Klinger 
1974]  whose  beam  is  reflected  from  a  mirror  attached  to  the  section  of  the 
cryostat  in  front  of  the  gradiometer. 

Figure  2.3  is  a  plot  of  the  differential  SQUID  response  at  twice  the 
shaker  drive  frequency  versus  the  autocollimator  tilt  signal  measured  at  the 
drive  frequency.  The  SQUID  and  autocollimator  responses  were  synchronously 
detected  with  a  PAR  5204  lock-in  amplifier.  The  upper  points  fit  the 
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expected  behavior  of  a  twice  frequency  signal  increasing  as  the  square  of 
the  tilt  amplitude.  This  plot  gives  a  scale  factor  of  1.2  +  .15  yV/E  with 
5.0  A  stored  in  the  differential  loop  and  a  drive  frequency  of  4.15  Hz. 

Equation  [l.33]  allows  us  to  predict  a  scale  factor  of  1.46  yV/E 
using  the  experimental  value  of  y  from  section  2.1,  and  a  baseline  of  3.2 
cm,  giving  reasonable  agreement  with  the  experimental  result. 

In  order  to  calibrate  the  system  with  a  gravitational  gradient  we  have 
measured  the  response  of  the  gradiometer  to  a  rotating  mass  quadrupole. 
Figure  2.4  shows  the  layout  of  this  experiment.  Two  masses  are  mounted 
on  an  aluminum  arm  a  distance  r  from  the  z  axis.  This  arm  rotates  in  a 
vertical  plane  about  the  z  axis  at  a  frequency  w.  The  gradiometer  is  posi¬ 
tioned  with  its  sensitive  axis  horizontal,  and  its  center  on  the  z  axis  a 
distance  R  from  the  center  of  the  mass  quadrupole.  The  azimuthal  angle 
is  the  angle  between  the  z  axis  and  the  gradiometer  axis. 

To  calculate  the  gravity  gradients,  note  that  the  displacement  vectors 
between  and  p^  and  M2  and  p^  are 

d^  =  (h sin  9  -r  cosojt)i  +  r  sinutj  -  (R  +h  cos0)  k  ,  [  2  4a] 


d2  =  (h  sinQ  +  r  cosoo  t)i  -  r  sinojtj  -  (R  +h  cos0)  k  . 


The  gravitational  acceleration  at  p  is 


[ 2.4b] 


To  find  the  gradient  we  project  this  onto  the  unit  vector  in  the  direction 


of  the  gradiometer  axis  and  take  the  derivative  with  respect  to  h.  So 
g(t)  =  ^a^t)  •  (cos9  k  -  sin0  i  )j. 


[  2.6] 


Taking  the  dot  product 


„w  d  f  -sin8  (h  sin  9+  r  coscot)+cos8  (R  +  h  cos9) 

o  v  t-  /  ,,  4  99?  9-13/2 
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[2.7] 


The  experiment  whose  description  follows  is  sensitive  only  to  time  varying 
gradients  so  in  taking  the  derivative  we  will  drop  terms  with  no  time  depen¬ 
dence.  This  leaves 
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2  2  2  2 
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Finally  we  set  h  =  0  and  again  drop  constant  terms.  This  leaves 


g(t) 


2  2 

3GMr  sin  9cos2k)t 
[r2  +  R2]5/2 


[2.9] 


We  find  that  the  only  time  dependent  term  is  at  twice  the  mass  quadru- 

pole  spin  frequency.  This  is  desirable  since  most  of  the  mechanical  noise 

produced  by  the  rotor  is  at  the  spin  frequency.  The  signal  also  has  a  clear 

2 

signature  in  that  rotating  the  gradiometer  about  the  vertical  gives  a  sin 
modulation  of  the  signal.  Unfortunately,  the  dependence  on  R  is  a  disadvan¬ 
tage.  For  R  >  r  the  signal  decreases  like  R  ^  so  the  rotor  has  to  be  rather 
close  to  the  instrument  and  relative  gradient  error  is  five  times  the  rela¬ 
tive  error  in  R.  For  these  reasons  the  rotating  mass  quadrupole  is  not  very 
useful  for  determining  the  scale  factor  of  the  instrument,  but  it  is  useful 
as  a  demonstration  that  the  gradiometer  is  actually  able  to  see  a  gravita¬ 
tional  gradient. 

Figure  2.5  is  a  photograph  of  the  mass  quadrupole.  The  masses  consist 
of  four  lead  bricks  mounted  with  their  centers  of  mass  35.6  cm  apart.  For 
this  experiment  the  gradiometer  dewar  was  hung  in  the  corner  of  an  acousti¬ 
cally  insulated  room  from  the  usual  latex  tubing.  The  mass  quadrupole  was 
positioned  just  on  the  other  side  of  the  wall.  With  this  arrangement  the 
dewar  could  be  rotated  to  change  9.  The  measured  values  of  the  parameters 
in  equation  [2.9]  were 


and 


M  =  43  kg, 

R  =  .58  m, 

r  =  .18  m, 

oj  =  20.1  rad /sec. 
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From  these  we  calculate 


3GMr‘ 


[R2  +  r2J5/2 


=  3.4  E  . 


The  mass  quadrupole  was  rotated  at  3.2  Hz,  and  the  signal  at  6.4  Hz  was 
synchronously  detected  with  a  lock-in  amplifier.  The  lock-in  reference  was 
provided  by  a  photodiode-phototransistor  switch  on  the  frame  of  the  mass 


quadrupole.  Figure  2.6  is  a  plot  of  the  gradiometer  signal  versus  azimuthal 

2„ 


angle  0.  The  solid  curve  is  a  fit  to  the  data  assuming  a  sin  9  dependence 
on  azimuthal  angle.  The  data  is  somewhat  noisy  because  electrical  noise 
from  the  mass  quadrupole  drive  system  causes  the  SQUID  to  unlock  frequently. 
When  the  SQUID  unlocks  its  output  changes  by  a  multiple  of  20  mV  causing  the 
lock-in  to  overload.  In  order  to  collect  data  a  short  averaging  time  had  to 
be  used.  The  time  constant  of  the  lock-in  amplifier  was  set  at  10  sec  for 
the  data  in  figure  2.6.  Using  the  data  from  the  previous  calibration  gives 


3GMr  ‘ 


U2  +  rV/2 


2.5  E  . 


This  is  significantly  less  than  the  calculated  value,  but  a  i  cm  error  in 

measuring  R  gives  an  error  of  about  10%  in  the  gradient.  In  addition  to  this 

error,  finite  size  effects  can  also  reduce  the  observed  signals.  The  fact 

2 

that  the  reponse  varies  like  sin  9  strongly  suggests  that  it  is  gravitational 
rather  than  seismic  or  acoustic.  We  also  verified  that  increasing  R  by  30 
cm  eliminated  the  signal  with  the  averaging  time  used  above. 

Figure  2.7  shows  the  gradiometer  output  noise  versus  frequency  between 
zero  and  ten  hertz.  The  large  peak  at  1.2  Hz  is  the  vertical  resonance  of 
the  suspension.  Above  2.5  Hz  the  noise  floor  is  at  3.3E//Hz,  although  there 
are  several  peaks  where  the  noise  is  more  than  40E//Hz\  Under  the  suspen¬ 
sion  resonance  the  noise  climbs  to  about  100E//Hz\  Two  intrinsic 
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Figure  2.6  Gradiometer  signal  from 
the  rotating  mass  quadrupole. 


sources  of  noise  are  readily  identified.  The  noise  performance  of  the 
commercial  SQUID  used  in  the  differential  readout  was  measured  separately 


before  we  installed  it  in  the  readout  circuit.  Between  zero  and  ten  hertz 
the  output  noise  was  1.1  pV//Hz  or  the  equivalent  of  .92  E//Hz  with  5.0  A 
in  the  readout  loop.  The  measured  performance  is  only  3.6  times  this  above 
2  Hz.  So  between  2  Hz  and  10  Hz  the  instrument's  performance  is  approaching 
one  of  the  intrinsic  limits  of  the  system. 

The  other  intrinsic  noise  source  is  the  Brownian  noise  in  the  differen¬ 
tial  mode  due  to  the  dissipation  in  the  system.  When  the  system  is  balanced 
the  differential  mode  is  just  a  harmonic  oscillator  subject  to  a  fluctuating 
force  with  power  spectral  density 

S  (to)  =  4KTH 

[2.10] 

according  to  Nyquist's  theorem.  Where  k  is  Boltzman's  constant,  T  is  the 

temperature,  and  H  is  the  mass  divided  by  the  energy  decay  time  of  the 

mode.  The  measured  energy  decay  time  of  the  differential  mode  with  5.0  A 

stored  was  6.3  seconds.  The  spectral  density  of  gradient  fluctuations 

corresponding  to  Sj(w)  is  given  by 

s  =  ✓Sf(u) 
g  -  , 

Mb  [2.11] 

at  frequencies  well  below  the  resonant  frequency  of  the  differential  mode. 

At  T  =  4.2°  K  this  gives 

S  =  0.18  E//Hz\ 


inn 


This  is  nearly  20  times  less  than  the  observed  noise  floor  of  the  instrument 
and  is  six  times  less  than  the  limit  of  the  SQUID. 

Although  the  noise  performance  above  2  Hz  is  quite  close  to  expectations 

the  noise  amplitude  between  0  and  1  Hz  remains  a  serious  problem.  In  order 

to  carry  out  the  inverse  square  law  test  we  need  to  measure  gradients  at  the 

center  of  a  cylindrical  shell  whose  central  horizontal  gradient  is  6  E. 

To  complete  a  highly  sensitive  test  we  would  like  to  measure  this  to  an 

_2 

accuracy  like  .5%.  This  requires  the  ability  to  detect  3  x  10  E,  a  signal 
well  below  the  noise  level.  However,  there  are  several  improvements  which 
may  reduce  the  noise  at  low  frequencies.  In  section  2.4  we  will  discuss 
the  sensitivity  of  the  instrument  to  thermal  fluctuations.  It  turns  out  that 
this  sensitivity  is  quite  high,  and  the  implementation  of  more  stable  thermal 
control  in  the  instrument  may  improve  the  low  frequency  noise  performance. 

We  will  also  complete  a  feedback  loop  to  stabilize  the  positions  of  the 
proof  masses,  and  this  may  reduce  low  frequency  noise  caused  by  nonlinearity 
in  the  suspension.  Such  nonlinearity  may  cause  rectification  of  high  fre¬ 
quency  signals ,  which  would  show  up  as  an  increased  noise  level  near  zero 
frequency. 

2 . 3  Feedback  Damping 

During  the  calibration  experiments  described  in  section  2.2,  feedback 
damping  was  used  to  greatly  reduce  the  signal  at  the  differential  mode 
frequency.  This  was  essential  because  excitation  of  the  differential  mode 
would  cause  the  SQUID  to  unlock  frequently  if  the  mode  was  left  undamped. 

Figure  2.8  shows  the  circuit  used  to  condition  the  feedback  signal. 

It  consists  of  a  bandpass  filter  with  adjustable  Q  and  center  frequency, 
a  phase  shifter  to  provide  90  degrees  of  phase  shift  at  the  bandpass  center 
frequency,  and  a  voltage  to  current  converter  which  provides  the  feedback 
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current.  According  to  equation  [1.62]  the  minimum  Qm^n  is  approximately 


1  -  co^/ov 


A 1  • 


[  2.12] 


to  can  be  measured  directly  since  it  is  the  normal  mode  frequency  in  the 
absence  of  damping.  We  can  calculate  an  expected  value  for  to^  from  equa¬ 
tion  [1.57].  According  to  [1.56]  and  [1.57] 

2  2  D (to  2  -  to  2) 

to„  =  to  +  1  o 

2  o  - - -  • 

(L1°  +  L2)(L3  +  L4)  [2.13] 


With  6.0A  in  the  readout  loop  we  had 


to^  =  505  rad/sec. 


Using  the  design  values  for  the  other  parameters 


D  2  2700  (pHr  , 

(L  °  +L2)(L3  +  L4)  =  4500 (pH) 2  , 


to  =  377  rad/s  , 

o 


we  find  an  expected  value  for  to2  of 
to2  =  458  rad/s. 

According  to  [2.12]  this  gives  Q  .  =11.  Using  [1.63]  we  find  the  Q  . 

min  mm 

should  be  reached  at 


WL1  +L2>  .  !  , 

-  =  to.  +  —  (to  -co„)  , 

n  12  1  /  ’ 


[  2.14] 


R  _  ,  =  3,1  x  10  ft  . 


According  to  [1.691 


M,  G  G  Q,  CD, 
f  v  s  f  f 


where  =  10  ^  H 

G  =  1.72  x  105  V/A 

s 

Qf  =  50 

CDj  =  500  rad/sec. 

With  these  values  of  the  parameters 

R  ,  ,  =  G  (4.1  x  103  V  -  fi/A) . 

ef  f  v 


[2.15] 


[  2.16] 


Figure  2.9  shows  the  measured  Q  of  the  differential  mode  versus  G  , 

v 

the  current  to  voltage  scale  factor  of  the  final  stage  of  the  feedback 


circuit.  As  G^  is  changed,  (d^  must  be  readjusted  so  that  the  bandpass 

filter  is  tuned  to  the  observed  frequency  of  the  mode.  If  Gy  is  made  very 

large,  the  frequency  falls  to  73.8  Hz  or  464  rad/sec.  This  is  quite  close 

3 

to  the  expected  value  of  cd^.  The  initial  Q  of  the  mode  was  3.2  x  10  . 

The  minimum  Q  obtained  was  close  to  50  rather  than  11.  This  is  probably  due 
to  the  Qj  of  the  bandpass  filter  in  the  feedback  circuit  which  was  left  at 
50  as  G^  was  varied.  The  value  of  Gv  which  gives  a  minimum  Q  is  larger 
than  expected  from  [2.14]  and  [2.16] ,  but  this  could  easily  be  due  to  a 


smaller  value  of  m^  than  the  design  value  of  1  pH. 

In  operation  this  feedback  loop  reduced  the  signal  at  the  resonant 
frequency  by  a  factor  of  approximately  30  and  greatly  improved  the  stability 
of  the  SQUID. 


2.4  Thermal  Sensitivity 


The  readout  described  in  section  1.2  is  sensitive  to  changes  in  the 


temperature  of  the  gradiometer  due  to  the  variation  of  the  penetration 
depth  in  niobium  with  temperature.  A  magnetic  field  applied  to  the  surface 


of  a  superconductor  will  decay  exponentially  inside  the  superconductor 
in  a  characteristic  length  X,  the  penetration  depth.  This  penetration 
depth  is  a  function  of  temperature,  and  it  is  part  of  the  gap  between  the 
readout  coil  and  the  opposite  proof  mass.  So,  as  it  changes  with  tempera¬ 
ture  the  readout  gap  varies  and  a  signal  is  detected.  An  approximation 
to  the  temperature  dependence  of  X  is  given  by  the  two-fluid  approximation 


t  2.17] 


[  Tinkham  1975]  ,  where  X ^  =  440Aand  t^  =  9.2°  K  for  niobium.  Taking  the 

derivative  of  \(t)  with  respect  to  temperature  gives 

3r  4,  -3/2 


dA 


[ 2.18] 


At  4.2  K 


=  9.7  A/°K  . 

d  t 


[2.19] 


Looking  back  at  the  introduction  we  see  that  a  gravity  gradient  of 
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1  E  produces  a  Ax  of  about  2  x  10  cm,  or  2  x  10  A,  so  we  can  expect 
a  small  change  in  temperature  to  produce  a  large  signal.  This  effect  is 
larger  than  the  length  change  caused  by  the  coefficient  of  thermal  expan¬ 
sion  of  the  niobium.  At  20°  K  this  coefficient  is 


1 

l 


dL 

dt 


3  x  10"7/°K 
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[Corruccini  1961] ,  and  it  is  at  least  one  order  of  magnitude  smaller  at 
4.2°  K.  Using  the  size  of  the  gap  between  the  proof  masses  for  L  we  have 

^  <  .4  A/°K  . 

dt 

Measurements  of  the  thermal  sensitivity  of  the  instrument  as  a  function 
of  stored  current  in  the  readout  coil  and  gradiometer  temperature  have  been 
carried  out.  Figure  2.10  shows  dV/dt  at  4.5°  K  versus  readout  current, 
where  V  is  the  output  voltage  of  the  differential  SQUID.  dV/dt  increases 
linearly  with  Iq,  and  the  slope  of  the  fit  is 

31.6  V/°K-A. 

The  measurements  were  made  with  a  thermometer  and  a  temperature  controller 
which  stabilized  the  temperature  by  heating  a  coil  wrapped  around  the  case 
of  the  gradiometer.  The  thermometer  was  a  carbon-glass  resistor  model 
CGR-1  manufactured  by  Lake  Shore  Cryotronics.  Its  conductance  was 
measured  with  a  potentiometric  conductance  bridge  made  by  S.H.E.  Corporation, 
and  the  temperature  was  controlled  with  a  temperature  controller  also  made 
by  S.H.E.  The  response  of  the  SQUID  to  temperature  changes  of  a  few  milli- 
degrees  was  measured  at  each  current,  and  no  hysteresis  was  apparent  in 
the  response.  The  slope  above  is  about  9  times  larger  than  would  be  pre¬ 
dicted  by  considering  dA/dt.  With  5A  stored  the  scale  factor  is  1.2  UV/E 

-14 

and  IE  produces  a  relative  motion  of  1.4  x  10  cm,  so  the  scale  factor 
for  displacement  is 

8 . 6  x  10^  V/cm. 

Using  [2.19],  and  assuming  that  the  magnetic  field  penetrates  the  readout 
coil  and  the  proof  mass  equally  we  get 

8.1  x  10^  V/cm. 
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Figure  2.11  Temperature  sensitivity  as  a 
function  of  temperature. 
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Part  of  this  discrepancy  could  easily  be  due  to  penetration  depth  changes 
in  other  parts  of  the  circuit.  For  instance  the  inductance  of  the  trans¬ 
former  primary  will  change  with  temperature.  It  is  also  possible  that  the 
effective  area  of  the  wire  in  the  readout  coil  is  larger  than  the  area  of 
the  proof  mass,  causing  a  larger  change  in  inductance  with  temperature. 

Figure  2.11  shows  the  behavior  of  the  dV/dt  as  a  function  of  tempera¬ 
ture.  The  measurement  points  were  taken  with  100  mA  in  the  readout  loop. 

The  solid  curve  shows  the  behavior  of  dA/dt  as  a  function  of  temperature 
from  [2.18]  .  This  curve  is  normalized  to  cross  the  measurement  at  4.7°  K. 
The  measured  points  do  not  increase  as  rapidly  as  [ 2.18]  but  the  rapid 
increase  in  dV/dt  as  temperature  approaches  9.2°  K  does  indicate  that  the 
temperature  dependence  is  related  to  the  superconductivity  of  niobium. 

This  very  large  temperature  sensitivity  may  explain  the  excess  gradio- 
meter  noise  below  2  Hz.  With  5A  stored,  the  thermal  sensitivity  is  158  V/°K 

compared  with  a  scale  factor  of  1.2  pV/E,  so  temperature  fluctuations  must 
—8 

be  kept  below  10  °K  to  prevent  them  from  contributing  to  the  observed 
noise.  We  currently  have  a  temperature  control  loop  with  a  stability  around 
10  k  °K,  and  the  excess  low  frequency  noise  is  about  one  hundred  times  the 
SQUID  noise. 

2 . 5  Flux  Creep 

Early  experiments  with  superconducting  gradiometers  [  Paik  1978]  indi¬ 
cated  that  readouts  of  the  type  described  here  would  show  drift  caused  by 
flux  creep  if  type  II  superconductors  were  used  in  the  readout  coils.  This 
drift  is  caused  by  the  motion  of  trapped  flux  lines  in  the  magnetic  field 
gradient  of  the  readout  coil  [Anderson  1962]  ,  [Anderson  1964]  .  Flux  is 

trapped  in  type  II  superconductors  when  H  ,  of  the  material  is  exceeded. 

cl 


The  early  readouts  were  constructed  using  Nb-Ti  wire  which  has  a  very  low 

H^.  These  readouts  showed  a  steady  drift  like  that  expected  from  flux 

creep  [  Kim  1962] .  To  eliminate  this  problem  we  investigated  the  stability 

of  pure  niobium  wire  which  has  H  ,  -  1000  G. 

cl 

In  order  to  reduce  the  time  required  for  cooldowns  and  changes  in  the 
experiment  we  constructed  a  small  probe  to  operate  as  a  dip  stick  in  a 
narrow  neck  storage  dewar.  The  probe  is  equipped  with  a  small  solenoid 
of  the  wire  sample  to  be  tested,  a  hall  probe  to  detect  the  field  produced 
by  the  persistent  current  in  the  sample,  and  a  pickup  coil  which  is  wrapped 
around  the  wire  sample  and  coupled  to  a  SQUID  to  detect  small  changes  in 
the  persistent  current. 

Figure  2.12  shows  the  drift  in  the  SQUID  output  with  2.0  A  stored  in 
a  sample  of  .\Tb-Ti  wire.  Figure  2.13  shows  the  drift  caused  by  5.0A  stored 
in  a  sample  of  pure  niobium  wire.  The  drift  in  the  niobium-titanium  is 
logarithmic  as  predicted  for  flux  creep  by  the  Anderson-Kim  theory.  The 
niobium  sample  does  not  show  any  steady  drift,  and  the  stability  of  the 
current  is  at  least  25  times  better  than  it  is  in  the  niobium-titanium. 
Because  of  these  results,  all  of  the  coils  in  the  gradiometer  were  built 
with  pure  niobium  wire  of  the  tvpe  tested  in  figure  2.13. 

2 . 7  Magnetic  Sensitivity 

Since  the  source  for  the  inverse  square  law  test  is  a  steel  cylinder, 
careful  magnetic  shielding  will  be  required  to  prevent  the  magnetic  field 
of  the  source  from  causing  a  spurious  signal.  Because  of  this  we  have 
built  several  layers  of  magnetic  shielding  into  the  cryostat.  The  instrument 
itself  is  the  first  layer  since  it  is  a  solid  superconductor  which  shields 
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the  pickup  coil,  transformers,  and  signal  leads.  The  inside  of  the  vacuum 
can  is  lined  with  another  lead  shield,  and  there  is  another  lead  shield 
and  a  j_i -metal  shield  outside  the  vacuum  can. 

This  system  was  tested  by  applying  a  magnetic  field  to  the  outside  of 
the  dewar  with  a  four  foot  coil  and  an  audio  amplifier.  Fields  up  to  2  milli 
gauss  were  applied  and  measured  with  a  flux  gate  magnetometer  without  any 
reponse  from  the  differential  or  common  mode  SQUID.  The  system  can  easily 


be  tested  to  higher  fields  if  necessary  with  a  more  powerful  amplifier  to 


drive  the  coil. 


CHAPTER  3 


Evidence  for  the  Inverse  Square  Law 


3.1  Celestial  Mechanics 


The  original  and  strongest  evidence  in  favor  of  the  inverse  square 
law  of  gravitation  comes  from  observations  of  the  mechanics  of  motion  in 
the  solar  system.  The  most  stringent  limits  on  variations  in  G  come  from 
the  agreement  between  theory  and  observation  of  the  precession  of  the 
perihelia  of  several  bodies  in  the  solar  system  and  the  binary  pulsar 
(PSR  1913  +  16).  Additional,  but  weaker,  constraints  on  possible  varia¬ 
tions  in  G  can  be  derived  by  comparing  determinations  of  planetary  masses 
made  by  observing  the  planet's  effect  on  the  orbits  of  other  bodies  at 
different  distances  from  the  planet  [  Mikkelson  1977]. 

The  orbits  of  bound  particles  in  an  inverse  square  force  field  are 
closed  ellipses,  and  small  deviations  from  this  force  law  cause  orbital 
precession  such  as  that  caused  by  corrections  due  to  general  relativity. 
Consequently,  any  variation  of  G  with  distance  would  cause  a  discrepancy 
between  observed  rates  of  precession  and  those  predicted  from  the  interaction 
of  the  planets  and  relativistic  effects.  Consider  the  motion  of  a  body  of 
mass  m  in  a  central  force  field  F(r),  since  angular  momentum  is  conserved 


0  =  L 


[3.1] 

is  a  constant,  and  the  radial  equation  of  motion  becomes  [Goldstein  1950] 


r  -  L^/mr^  =  F(r)/m 


[3.2] 


For  a  nearly  circular  orbit  of  radius  a  we  can  expand  [  3.2]  in  a  power 
series  in  x  =  r  -a,  so 
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Using  [3.2]  to  cancel  the  terms  which  are  constant  we  have 


x  -  +  F'(a)j  x  =  0  , 


[3.4] 


to  first  order  in  x.  From  this  we  can  see  that  the  body  oscillates  harmoni¬ 


cally  about  the  circular  orbit  with  a  period 


T  =  2tt 


- [ 3F (a)+aF 


-  P 

'  ( a )]  j  ‘ 
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Since  3  is  approximately  constant,  the  advance  of  a  planet's  perihelion  in 


one  orbit  is 


19  =  T  9 


-  [ 3F (a)+aF  ’  (a)] 


Using  [3.2]  again 


Ad  =  2tt 
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Mow  suppose  that  G  does  vary  with  r,  so  that 


Applying  [3.7]  we  find 
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'  [ 3.8] 
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to  first  order  in  G'(a)/G(a).  The  most  accurate  measurements  of  orbital 
precession  have  been  done  on  the  orbit  of  Mercury  where  the  general  relativ¬ 
istic  effects  are  largest.  The  predicted  value  of  A9  for  Mercury  after 
subtracting  out  the  effects  of  other  planets  on  its  orbit  is 

A0  =  5.01  x  10  7  rad/orbit  , 
mer 

where  A0  is  calculated  with  the  relativistic  formula  [Marion  1970] 

A9  =  6G77  Mg/[c2a(l  -  e2)3  , 

where  Mg  is  the  solar  mass,  C  is  the  speed  of  light,  a  is  the  semimajor 

axis  of  the  orbit,  and  e  is  the  orbital  eccentricity.  The  measured  value  of 

A9  is  [  Shapiro  1976] 
mer 

A0  =  5.03  +  .02  x  10“7  rad/orbit, 

mer  — 

The  excellent  agreement  of  these  numbers  places  a  stringent  limit  on 
aG'(a)/G(a)  for  a  =  5.8  x  107  km,  the  semiraajor  axis  of  Mercury's  orbit. 
Measurements  on  the  orbit  of  Mars  [  Laubscher  1971]  and  the  orbit  of  Icarus 
[Shapiro  1971]  also  give  good  agreement  with  general  relativity.  These 

g 

measurements  rule  out  variations  of  G  of  more  than  a  few  parts  in  10  at 

g 

ranges  close  to  10  km.  Observations  of  the  binary  pulsar  PSR  1913  +  16 

[Taylor  1976]  indicate  that  G  is  constant  to  parts  in  10^  at  10^  km. 
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A  comparison  of  G  at  distances  down  to  about  10  km  with  G  at  10  km 
can  be  made  by  comparing  determinations  of  the  masses  of  the  inner  three 
planets  made  at  these  two  distances.  The  first  mass  determination  comes 
from  analysis  of  inner  planetary  radar  ranging  data  where  the  distances 
between  the  planet  and  the  body  whose  orbit  it  perturbs  is  on  the  order 


of  10  km  (Shapiro  1973]  .  For  Mercury  and  Venus,  the  second  mass  determina¬ 
tion  comes  from  observations  of  the  trajectory  of  Mariner  10.  Mariner  10 
flew  by  Venus  with  a  minimum  separation  of  10^  km  [Howard  1974]  and  came 
within  3  x  10  km  of  Mercury  [  Esposito  1976]  .  The  second  mass  determination 
for  the  earth-moon  system  was  made  with  lunar  laser  ranging  data.  This 

3  o 

comparison  leads  to  a  limit  on  the  change  in  G  between  3  x  10  km  and  3  x  10 
km  of  5  x  10  4  [  Mikkelson  1977] . 

In  summary,  the  inverse  square  law  is  very  well  verified  at  astronomical 

distances  and  any  proposed  violation  must  be  extremely  small  at  these  ranges. 

However,  the  exponential  deviation  of  equation  [ 1 . 1]  has  this  property  and 
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is  not  constrained  by  solar  system  data  if  $  <<  10  km.  Long  has  proposed 

a  violation  of  the  form 


G(r)  =  G(1  +  .002  2, nr) , 


[3.10] 


[Long  1976]  based  on  laboratory  data,  but  some  mechanism  must  cut  such  a 
violation  off  at  distances  of  1000  km  because  it  predicts  a  change  in  G  of 
10%  between  10^  km  and  10^  km. 

3 . 2  G  Measurements 

At  laboratory  distance  scales  measurements  of  G  provide  some  slight 
evidence  that  G  is  constant.  Long  has  reviewed  a  number  of  precise  labora¬ 
tory  measurements  [  Long  1974]  .  A  plot  of  these  measurements  is  shown  in 
figure  3.1.  The  measurements  are  consistent  at  the  level  of  three  standard 


deviations,  but  Long  points  out  that  variations  in  G  of  less  than  1%  over 
this  distance  range  are  allowed  by  the  data.  He  fits  the  dashed  line  to 
the  data  using  equation  [3.10],  and  his  own  measurements  at  4.5  cm  and 


Most  of  the  measurements  in  figure  3.1  were  made  with  torsion  balances; 
the  technique  pioneered  by  Cavendish  in  1798,  although  Poynting  and  Richarz 
used  ordinary  balances  to  detect  the  change  in  weight  on  a  test  mass  from 
an  attracting  mass.  The  most  modern  measurement  of  Rose  et  al.  was  made 
by  measuring  the  angular  acceleration  on  a  torsion  balance  mounted  on  a 
rotating  table  subject  to  a  constant  torque  from  masses  also  mounted  on  the 
table.  The  references  for  the  various  points  are: 

B  [  Boys  1895] 

BR  [  Braun  1896] 

P  [  Poynting  1891] 

RK-M  [Richarz  1898] 

H  [Heyl  1930] 

HC  [Heyl  1942] 

R  [Rose  1969],  also  see  [Beams  1971]. 

Long's  measurement  also  employed  a  torsion  balance.  He  used  two  rings 
as  the  attracting  masses  and  positioned  the  test  mass  of  the  torsion  balance 
at  the  point  where  the  derivative  of  force  with  respect  to  distance  is  zero 
for  each  ring  in  order  to  reduce  the  necessary  positioning  accuracy.  Although 
the  G  measurements  cannot  be  construed  as  verifying  Long's  measurement,  they 
do  leave  open  the  possibility  that  something  very  fundamental  is  going  on. 

Long  has  pointed  out  another  set  of  measurements  which  are  consistent  with 
his  own.  These  were  made  by  measuring  the  frequency  shift  in  a  torsion 
balance  as  an  attracting  mass  was  moved  between  several  locations  [  Karagioz 
1976] .  Long  has  analyzed  their  data  and  found  a  value  for  X  in  the  equation 


G(r)  =  G [l  +  X  In  r]  , 


X  =  ^2.01  +  .6l)  x  10' 


[  3.11] 


[  Long  1979]  . 


3.3  Recent  Experiments 


Since  1976,  when  Long  first  claimed  that  a  violation  of  the  inverse 
square  law  exists,  some  work  has  been  done  with  the  specific  purpose  of 
checking  his  results.  The  only  work  to  claim  an  accuracy  better  than  Long's 
was  done  at  U.C.  Irvine  [  Spero  1980]  using  a  method  similar  to  that  described 
here.  The  U.C.  Irvine  group  used  a  torsion  balance  to  measure  the  horizontal 
gravitational  gradient  at  the  center  of  a  stainless  steel  cylinder.  Their 
cylinder  has  a  diameter  to  length  ratio  of  about  .1  and  is  60  cm  long.  The 
gradient  was  measured  by  moving  the  cylinder  so  that  the  proof  mass  of  the 
torsion  balance  was  next  to  one  wall  and  comparing  this  to  the  force  on  the 
balance  with  the  proof  mass  next  to  the  other  wall.  They  used  y-ray  trans¬ 
mission  studies  to  measure  the  homogeneity  of  their  source.  They  find  a 
value  for  X  in  [3.11]  of 


V  * 
•V- kV.S 


•  .  J 


i.  .  J 


P  .  .  ) 


PL  .  i 


.  J 


X  =  (1  +  7)  x  10  . 

This  would  appear  at  least  to  rule  out  Long's  logarithmic  model  for  the  viola- 


t 


tion,  but  Long  has  argued  that  such  a  form  might  arise  from  a  vacuum  polar¬ 
ization  effect  similar  to  the  effect  which  produces  a  logarithmic  modifica¬ 
tion  of  Coulomb's  law  at  short  distances  [Long  1980a].  Such  an  effect  might 


» 


be  unobservable  in  a  null  experiment  due  to  the  lack  of  a  polarizing  field 


[Long  1980b] 
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Several  other  interesting  experiments  have  been  done,  but  none  has 
sufficient  accuracy  to  rule  out  Long's  result.  Panov  and  Frontov  [Panov 
1979]  used  a  torsion  balance  to  measure  the  ratio  of  G  at  3  m  to  G  at  .4  m, 
and  the  ratio  of  G  at  10  m  to  G  at  .4  m.  Their  results  were 


G(3m) 
G( . 4m) 


1.003  +  .006  , 


and 


G(10m) 
G( .  4  m) 


.998  +  .012  . 


Unfortunately,  the  stated  errors  are  comparable  to  the  size  of  the  effect 
predicted  by  Long.  Hirakawa  [  Hirakawa  1980]  has  measured  the  interaction 
between  a  gravity  wave  antenna  resonant  at  60.5  Hz  and  a  steel  bar  rotating 
at  30.3.  His  analysis  of  the  data  confirms  Newton's  law  to  +  3%  over  the 
range  from  2  to  5m,  but  this  error  is  almost  ten  times  the  expected  effect 
from  Long's  work.  Taken  as  a  whole,  this  work  indicates  that  more  precise 
measurements  at  different  ranges  will  continue  to  be  interesting  and 
challenging  for  some  time. 


1 


I 


t 


< 


In  this  chapter  we  will  describe  the  source  and  some  possible  signals 
from  our  planned  experiment.  Our  source  is  a  large  steel  cylindrical  shell 
which  is  designed  to  fit  over  the  dewar  containing  the  gravity  gradiometer. 


This  cylinder  was  machined  from  a  centrifugal  casting  of  ASTM-27  grade 

60-30  normalized  steel.  Its  length  is  5.080  m,  the  inside  diameter  is 

0.4572  m,  and  the  outside  diameter  is  .5588  m.  This  gives  it  a  mass  3.24 
3  3 

x  10  kg  or  7.15  x  10  lbs.  For  the  moment  we  will  assume  that  this  source 
is  perfectly  cylindrical  and  homogeneous  for  the  purpose  of  calculating  its 
gravitational  field.  We  will  analyze  the  errors  due  to  deformation  and 
innomogeneity  of  this  source  in  a  later  section. 

The  acceleration  produced  by  the  gravitational  field  on  a  test 
particle  positioned  at  p(r,<]>,z)  is  given  by 

L/  2  2T r  R+t 


a  (r  ,<J> ,  z) 


■/// 

-L/2  0  R-t 


—  G(^-~P')  r'dr’dtjj'dz’ 

\t-n3 


[  4.1] 


where  G  is  the  gravitational  coupling  constant,  a  is  the  density  of  the 
#  ■+ 

cylinder,  p’  is  the  integration  variable,  and  R  is  the  mean  diameter  of  the 
cylinder.  The  coordinates  are  shown  in  Figure  4.1.  It  is  clear  from 
symmetry  considerations  that  the  result  for  a  must  be  independent  of  $  so 
we  may  set  $ =  0  and  consider  only  a  test  particle  in  the  z-x  plane.  Then 


and 


_►  A  A  /V 

p  =  rr  +  zz  =  ri  +  zk  , 


A  A  A 

p '  =  r '[  i  cos -I)'  +  j  sin  4> ']  +  z '  k. 


So 

A  A 

p  -  p '  =  (r  -  r '  cos  4> ' )  i-  r  '  sin  <j> '  j  +  (z  -  2  *  )k  , 

[4.2] 

and  for  a  we  have 

L/2  2tt  R  +  t 

a(r,z)  =  -Go  f  f  f  [(r-r  *cos<j)  * )  i-rsind/j  +  (2  -2  '  )k]  r  'dr  1  d<j)  'dz 

J  I  J  7  2  0  3/2 

-L/2  0  R-t  [(2-z')  +r  +  r’  -2rr ’cos<i> ']  [43, 

Note  that  the  integral  over  (p  1  in  the  j  term  vanishes,  since  this  term  has 
odd  parity. 

Let  us  define  the  dimensionless  variables  2z/L  =  £,  2z’/L  =  r/R  =  5, 
r'/R  *  a,  f|  =  2R/L ,  and  a)  =  t/R.  Then 

1  2TT  1+0} 

a(r,z)  =  -  Go  p2  L  ffl  a  [  ILiiz  a  cos(t>  '  )*  +  —  C  ' )  k]  dad^'dC'. 

2  1//  3/2 

-1  0  1-0)  [(C-C')2+n2(C?+  a2  -  2^acos<J)')] 

[4.4] 

Our  experiment  will  entail  measurements  of  gradients  of  the  cylinder, 
so  we  want  to  calculate  3az/3z,  and  3ar/3r.  We  also  need  to  know  how  these 
quantities  vary  with  position  near  the  center  of  the  cylinder,  in  order  to 
determine  how  accurately  the  gradiometer  needs  to  be  positioned  in  the 


cylinder.  First  we  calculate  3a  /3r  for  z  =0.  Then 


1  2tt  l+o) 


ar(r,o) 


-Gan  2l 


fff  a  n  (C  -  qcostp  *  )dad<}>  'dc  ' _ 

J  J  J  f  2  2  2  7  "1 

-1  0  1-co  k'  +n  (5  +  a  -  2£ctcos<j> ’ )J 


Doing  the  integral  over  we  have 


2ir  l+o) 


a^fr.o)  =  -2GOR 


//  q(C~  acos4> 1 )  [l  +  n2  (5  2  +  a2-2Cotcos(i) 1 )]  d<()'dq 

0  1-ti)  [«2  +  “2-«“c“*’]  f4 


Now  since  h  =0.01  for  our  source  we  can  expand  in  a  power  series  in  n 
and  then  integrate  over  term  by  term.  Then 


2tt  l+o) 


a^Cr.o)  =  -2GcR 


// 

0  l-o) 


o(g  -qcostp'-)  r 

[  £2+  q2  -2£qcosc(> '  ]  L 


_  n  ?  ?  3422  2 

2~  (£z  +  qz  -  2C3icos<J) ' )  +  g  n  (K  +  ct  -2£qcos<p') 


5  6  .,2  2 

16  n  +* 


5  2  +  a2  -2Cctcos<t')3  +  0(n8)]  . 


The  ■+>'  integration  over  the  first  term  gives  zero,  so 

1+0) 

a  (r,o)  =  2  G  OR  f  qn^cfl  -  3/4  n2  (£2  +  2q2) 

r  1-0)  l 

+  |  n4  [(C2+a2)2+2(C2+a2)q2+2qV]  -  ff  n6  [(C2  +  < 


+  3q2(C2  +  q2)  +  6  q2£2(£2  +  q2)  +  3  q4£2  I  dq  . 


The  cylinder  mass  is  given  by 


M  =  tt[(R+  t)2  -  (R  -  t)2"lLo  =  2R2Lttct  f 


Let  us  define 


/*  n  (I+03)  -  (1-0)) 

A  (03)  =  /  a  d  a  =  - - - 

n  I  n  +  1 


[4.10] 


Then  M  =  2R2LTOA1(o3>  .  Finally  we  do  the  integration  over  a.  We  have 


ar(r ,0) 


2  f  -,2  /.  2  A_  (03 )  \ 

GMn  ,  L  _  3n_  fe  +  _3 - ) 

IT  c  l  2  \  2  +  AL  (03)/ 


fr4  0  A  (03)  a  (03) 

15  4  K  ,  o  r2  _J3 _  ,  _ 

8  n  l  3  2  5  Ax(03)  A1(o3) 


35  6  [sl  .3  £Vl+i  52  Vl+  Vl|  +o(n8)i 

16  |  4  A^  (03)  2  A^(03)  A^(03)J  J 


[4.11] 


A3(o>) 

A^  (oj) 


=  1+03 


A  5  (o) ) 
A3  (oj) 


10  2  4 

1  +  —  03  +03 


A?  (03)  2  4  6 

— —  =  1  +  7CO  +70)  +  o)  , 

A^  (03) 


[4.12] 

Since  n  is  1/100  and  u)  is  1/20,  only  the  first  two  terms  in  this  expression 
give  measureable  contributions.  The  solution  has  two  important  advantages. 
The  lowest  order  contribution  to  Oa^/Or  is  GM/R  L(2R/L)  .  This  is  quite 
small  since  2R/L  is  small;  consequently  any  new  non-newtonian  force  makes 
a  proportionally  larger  contribution  to  the  central  gradient.  Ideally  the 
newtonian  term  would  be  zero,  and  we  would  be  doing  a  true  null  experiment. 
However,  this  would  require  an  infinitely  long  cylinder,  so  we  must  settle 
for  a  cylinder  which  attenuates  the  newtonian  gradient  by  a  factor  of  100. 
The  second  advantage  is  that  the  gradient  in  the  central  plane  of  the 
cylinder  is  nearly  independent  of  position. 
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Consider  the  lowest  order  term  in  the  force  equation  [4.12],  which 

is  of  the  form  a  =  Crr  where  C  is  a  constant.  In  the  x-y  coordinates  this 

has  the  form  a  =  c  (xx  +  yy) .  Hence  8a  /8x  =  8a  / 8y  =  c  anywhere  in  the 

x  y 

2  2 

central  plane.  The  lowest  order  correction  to  this  is  -3/4  ri  £  which 
is  very  small.  If  we  want  to  measure  the  central  gradient  to  0.1%  we  need 
only  have  C2<(.1)  4/3,  or  r<(.37)R  =  18.8  cm.  This  is  a  very  weak 


centering  requirement. 

Next  we  consider  variations  of  8a^/8r  with  z.  Using  equation  [4.4] 


we  have 


»r<*>  ■  tfU 


-Gcm2L 


1  2tt  l-o) 

Iff  4 

-1  0  l+o)  > 


'  )2+(an)2] 


-  2  2  2,, 
3u  a  cos 

[(?  -4’)2+(an)2 


dad4>'d4'  • 


[4.13] 


Doing  the  '  integral 

,  .  -GotTTL  n2 
gr(z)  =  — 


I  1+0) 


If  r(  7 — 7 - TP 

-1  1-0)  '  ks-c,')  +(an) 


f(c  -C')2+(an)2 


dad^  ’ 


[4.14] 


This  function  is  shown  in  figure  4.2.  It  is  evidently  quite  flat 
near  the  center  of  the  cylinder  and  falls  off  rapidly  outside  the  end. 

To  determine  how  accurately  the  gradiometer  needs  to  be  positioned  verti¬ 
cally  to  measure  the  central  gradient  we  expand  9a  /9r|  in  £  for  small 

r  lr=0 

£  using  [4.16]  .  This  gives 


GM  2  f.  ,, 
—  n  J 1  +  3 

LR  [ 


a  )  + 


2  2 

oftV)  + 


a(n2)j 


[4.18] 


So  the  gradient  error  produced  by  a  vertical  positioning  error  is  second 
order  in  the  vertical  displacement  of  the  gradiometer  from  the  center  of 
the  cylinder.  Symmetry  requires  that  the  error  be  some  even  power  of  £. 

3 

An  error  of  less  than  1  part  in  10  requires  that  the  gradiometer  be  within 
4.6  cm  of  the  center.  This  is  more  severe  than  the  radial  positioning 
requirement  but  still  is  easily  obtainable. 

There  is  one  further  positioning  error  which  needs  consideration. 
Newtonian  gravitation  obeys  Poisson's  equation  so  that  if  4>  is  the  gravi¬ 
tational  potential,  then  in  vacuum 


^24>  = 


9 jp  B2<p 

2  +  2  +  2 
9x  9y  9z 


0  . 


[4.19] 


For  a  perfect  cylinder,  symmetry  requires  that 


This  means  that  as  the  gradiometer  is  tilted  away  from  horizontal  up  to 
vertical  its  output  will  change  by  a  factor  of  -2,  hence  a  degree  of 
leveling  accuracy  is  needed.  To  determine  this  we  need  to  evaluate  the 
gradient  signal  at  the  center  of  the  cylinder  as  a  function  of  angle 
between  the  gradiometer  sensitive  axis  and  the  horizontal  direction. 

If  we  project  the  acceleration  onto  this  direction,  we  have 


where 


n  =  i  cos  9  +  k  sin  a  . 


The  differential  of  a  is 

y 


da9  = 


3x 


Si 

3z  J  , 


and  the  directional  derivative  of  a  along  the  n  direction  is  then  given  by 

y 


3  a  R  Sa¬ 
ga  =  a —  cos  0  +  5 -  sin  0  , 

s9  3x  3  z 


[A. 22) 


see  for  example  [Buck  1965].  But 


=  ^7  Jax  cos  0  +  a^sinBJ 


3x  |  x 
3a 


3a 


r  cos  0  +  r — —  sinB  , 
3x  dX 


and 


^a0  3  fa  cos  0  +  a  sin0l 

3Z  "  3Z  l  X  z  * 


3a 


3a 


r — —  cos  9  +  5 — —  sin  8  . 
3z  dz 


Evaluating  these  quantities  at  the  center  of  the  cylinder,  we  have 


3a  3a 

_ x  _  _ z 

3z  3x 


=  0 


since  a  =  0  along  the  z  axis  and  a  i  0  in  the  central  plane  of  the 

X  2 

cylinder.  Using  equation  [4.22]  , 


89  = 


3a 

„  3a 

_ X 

cos  8  +  r — —  sin  0 

3x 

3  z 

3a 

_ X 

Icos2  0-2  sin2  0  I 

3x 

l  1 

3a 

_ X 

3x 

|l  -  3  sin2o|  . 

[4.23] 


These  results  can  also  be  obtained  by  considering  a  similarity  transforma¬ 
tion  on  the  tensor  T  which  rotates  the  spatial  coordinates  with  respect  to 
the  cylinder  axis  [Goldstein  1959]. 


In  order  to  have  3  sin2 0  <  10  3  we  need  0  <  1° .  This  seems  a  rather  severe 
constraint.  Fortunately  gg  at  0  = 0  is  a  maximum,  so  that  by  measuring  g 
at  three  values  of  0  which  we  separate  by  a  known  A0 ;  9a  /9x  =  g  can 

X  X 

readily  be  found  by  interpolation. 


To  summarize,  the  radial  gradient  at  the  center  of  the  cylinder  is 


-3/2 


GM 

LR" 


(f)2  *  (f 


[4.24] 


and  the  axial  gradient  is 


-2GM 

LR2 


2  ?1"3/2 

(?)[-(?)]  . 


[4.25] 


from  [4.21] .  These  gradients  have  the  important  characteristic  that  they 

change  slowly  as  the  gradiometer  is  moved  away  from  the  cylinder  center. 

Careful  attention  needs  to  be  given  to  the  orientation  of  the  gradiometer 

axis  with  respect  to  the  cylinder  axis,  however. 

The  orientation  error  as  well  as  errors  due  to  perturbations  of  the 

cylinder  discussed  in  sec.  4.3  can  be  eliminated  by  planning  the  experiment 

to  measure  g  ,  g  and  g  .  The  gradiometer  mount  described  in  Chapter  1 
X  g  z 

allows  the  instrument  to  be  supported  in  two  positions  which  orient  the 
gradiometer  axis  vertically  and  horizontally.  Two  positioning  holes  in  the 
mount  ensure  that  the  axis  tilts  by  90°  between  these  positions.  Equation 
[4.19]  says  that  the  sum  of  g^,  g^  and  g ^  is  zero  if  the  inverse  square 
law  is  correct.  This  remains  true  In  the  presence  of  distortions  of  the 
cylinder  and  orientation  errors. 

However,  we  will  see  in  sec.  4.4  that  the  finite  size  of  the  gradio¬ 


meter  puts  a  limit  on  the  accuracy  with  which  equation  [4.19]  can  be 


verified. 


4 . 2  Gradients  produced  by  a  Yukawa  Potential  and  the  Logarithmic  Potential 


As  discussed  previously  several  authors  have  proposed  a  violation  of 
the  inverse  square  law  which  has  the  form  of  an  exponentially  damped  force. 
In  this  section  we  will  consider  the  gradient  in  our  cylindrical  source 
which  such  a  force  law  would  produce.  This  force  results  from  a  potential 
of  the  form 


(r)  „  -oGMe 
Yy  r 


-Sr 


[4.26] 


The  two  parameters  a  and  0  determine  the  magnitude  and  range  of  the  potential 


Since  <p  is  exponentially  cut  off  for  r>>  1/S,  we  can  only  expect  mass  at  a 


,-l 


range  of  S  to  produce  significant  variations  from  the  inverse  square  law. 
The  force  from  [4.261  is 


■*  .  .  -otGM  ..  ,  „  .  -Sr  -*• 

ay(r)  =  ~3  (1  +Br)e  r  •  [A. 27] 

r 

Using  this  we  can  write  an  equation  for  the  force  on  a  test  particle  in  the 
field  of  a  cylindrical  shell, 

L/2  2  tt  R+t 

a(r  z)  =  -aGo  f  J  f  r  '  [ (r-r  ' cosft  ' )  i  +  (  z-z  ' )k  ]  (1+Sp)c  8pdz ' d4>  Mr  '  > 
-L/2  0  R-t  P 

[4.28] 

where 

p  =  |(z-z  ' )  r  ^  +r  '  ^  -  2rr  1  cos  <p  ']  ^  * 


0  s  /v  \ 

We  would  like  to  compare  the  gradients  y^  =  — ^ — 


r=0 


_  9 (a-z) 
and  yz  "  —37“ 


r  =0 


with  the  results  for  the  Newtonian  gradients,  so  we  first  calculate 


dr' dz’ 


yr(z)  =  outgo 


j  j  n 

-L/2  R-t  q  1 


£_  1 3r  (1+Bq)  (3  r  ‘-2)  -  2Bq 


[4.29] 

2  2  % 

where  q  =  t(z-z')  +  r'  ]  .To  lowest  order  in  t/R  we  can  assume  that  r' 


in  the  integrand  takes  the  constant  value  R,  then 

1 


yr(z) 


2tanGo 


R 


+  (eV-2)-^ 


-BRs 


n 


df  , 

[4.30] 


2  2  ^ 

where  s  =  [(£-  £')  +  n  ]  .  Using  M  =  (2t)(2TTR)L, 


.  .  aGM 
y  (z)  =  — j 

2LR 


BRs 

e  n  dC'  . 


'-1 


[  4.31] 


Figure  4.3  shows  yr(0)  calculated  numerically  as  a  function  of  the  BR. 

When  the  range  B  1  is  comparable  to  the  cylinder  radius  there  is  an  enhance 
ment  of  the  gradient  of  almost  50  compared  with  what  a  Newtonian  force  law 
with  the  same  coupling  constant  would  produce.  As  B  becomes  small  the 
exponential  force  law  looks  Newtonian  and  we  recover  the  result  of  section 
4.1.  As  B  ^  becomes  much  less  than  R  the  force  along  the  cylinder  axis 
decays  exponentially  to  zero. 

Figure  4.4  shows  the  variation  of  y  with  z  for  three  values  of  BR. 

As  expected  if  BR  is  small,  the  signal  looks  approximately  like  that  of 
figure  4.2,  the  Newtonian  limit.  As  BR  becomes  about  1  the  signal  is 
large  and  independent  of  position  as  long  as  the  gradiometer  is  inside 
of  the  cylinder.  As  the  instrument  leaves  the  cylinder  the  signal  falls 


rapidly  Co  zero  due  to  the  exponential  factor.  The  signal  for  BR  =  5  is 

similar,  but  the  signal  inside  the  cylinder  is  much  smaller,  since  the  walls 

are  now  beyond  the  range  of  the  force,  (3  ^ . 

If  the  gradiometer  is  oriented  along  the  cylinder  axis,  we  can  calculate 

the  gradient  y  (z)  in  a  manner  similar  to  that  used  to  obtain  equation  [4.29] 
z 

Using  [ 4 . 26] , 


az(r ,z) 


-aGcr 


L/2  2n  R+t 

iff1 


-L/2  0  R-t 


(z-z  ’ ) 
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dr  'dtp  'dz  1 . 


[4.32] 


Since  we  are  interested  here  in  forces  along  the  z  axis  we  can  set  r  =  0  and 
do  the  integral  over  <p ' .  To  lowest  order  in  t/R  we  can  also  take  r'~  R 
and  do  the  integration  over  r'.  Then 


a  (o,z) 
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[4.33] 


So  the  gradient  is 


Y,(z> 
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A  numerical  calculation  of  y^/O)  as  a  function  of  8R  is  shown  in  figure  4.5. 

As  in  figure  4.3  as  8R  becomes  large  y  (0)  goes  to  zero  exponentially.  In 

z 

this  case  however,  the  fall  off  starts  at  8R  *  .1.  This  reflects  the  fact 
that  n  =  .1.  As  soon  as  the  range  8  ^  is  small  enough  that  the  missing  mass 


at  the  ends  of  the  cylinder  is  invisible,  each  proof  mass  sees  exactly  the 


same  mass  configuration  and  there  is  no  gradient.  As  3  ^  becomes  small  we 

,  2 


again  recover  the  Newtonian  limit  which  from  section  4.1  is 


-2an 


GM 


(i+n2)3/2  r2l 


y  (0)  is  quite  different  from  y  (0)  in  that  there  is  no  large  enhancement 
z  r 

of  the  gradient  between  these  two  limits.  Consequently,  the  experiment  must 
be  done  with  the  gradiometer  horizontal  to  see  a  deviation  from  the  inverse 
square  law. 

Unlike  the  Newtonian  potential,  the  Yukawa  potential  does  not  obey 
Poisson's  equation,  and  there  is  no  simple  relationship  between  y^CO)  and 
y^CO)  from  all  values  of  3.  Suppose  that  a  potential  of  the  form  [4.26] 
does  exist  with  the  optimum  range  for  this  experiment.  This  is  the  case 
for  3R  —  1.5  or  S  2  -  16.7  cm.  From  figure  4.3 


,  .  ,  .  /Q  2  GM 

y  (o)  =  y  (o)  =  48  n  a  — =■  . 
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From  figure  4.5  y  (0)  is  nil  in  this  case  so  the  sum  of  the  gradients 


,  w  .  ,  ,  .  ,  2  GM  ~  aGM 

y  (°)  +  v  (o)  +  y  (o)  =  96  a  n  — ~  ~  — ~  . 

y  z  LR  LR 


4.36] 


An  alternative  to  the  Yukawa  potential  has  been  proposed  by  D.  Long 
[  Long  1974]  .  It  is  a  logaritlimic  potential  of  the  form 


r  \  GM  0 

V  (r )  =  -  a  —  in 


M 


4.37] 


where  r^  is  a  range  which  Long  takes  to  be  1  cm.  Using  this  potential 


-*  ,  Ln(r/r  ) 

a  (r )  =  -VV  (r)  =  a  GM  — - - j-2- 

r  r 


r  . 


[4.38] 
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The  first  term  has  exactly  the  form  of  the  Newtonian  force  scaled  by  a 

-3 

factor  of  a.  Since  a  is  presumably  small,  (Long  gives  cx  =  2  x  10  )  we 
ignore  this  term  in  what  follows.  As  before,  we  would  like  to  calculate 
the  gradients  which  result  in  a  cylinder.  So  we  write  the  acceleration  on 


a  test  particle  in  the  cylinder 
L/  2  2-n  R+t 


a(r  ,z)  =  -O.CO  J  f  J  r'[(r-r'cos<}>'i+  (z-z')k]Ln(p/ro)  dr,d(£,dzI> 
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For  comparison  with  figures  4.3  and  4.5,  let  us  first  calculate  L^(z)  = 

3a  /3rj  .  Taking  the  derivative  first 
r  '  r=0 
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Setting  r  =  0  and  doing  the  integration  over 
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to  lowest  order  in  R/t.  Making  the  usual  change  of  variable  <;'  =  2z'/L 


we  find 
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Similarly  we  can  calculate  L^fz)  to  be 

Lz(z)  =  ^tn2  ^2j|[£n(sR/nro)  +  (c-C')2/s2]  df.  [4>43] 

A  numerical  calculation  of  L^CO)  and  Lz(0)  is  shown  in  figure  4.6  for  a 

wide  range  of  R/rQ.  The  solid  curve  is  L^CO),  and  the  dashed  curve  is 

L  (0) .  As  for  the  Yukawa  potential,  L  (0)  is  much  greater  than  g  (0)  if 
z  r  r 

ct  —  1 ,  And  again  L^CO)  <<Lr(0).  In  this  case  however,  there  is  no  Newtonian 
limit  and  Lr(0)  is  only  very  weakly  dependent  on  the  choice  of  range  r  . 

Using  Long's  choice  of  a  and  r^  [Long  1976]  we  find  that  Lr(0)/gr(0)  =  -.185, 
a  very  large  effect.  This  result  has  effectively  been  ruled  out  by  the 
experiment  of  Newman  [Spero  1980]  .  However,  Long  has  recently  suggested  that 
a  logarithmic  modification  of  the  inverse  square  law  might  not  be  observable 
inside  of  a  cylinder  [ Long  1980b] . 

4 .3  Gradient  Errors  due  to  Cylinder  Perturbations 

Among  the  principal  sources  of  error  in  a  measurement  of  the  cylinder's 
gravity  gradients  are  perturbations  of  the  cylinder's  density  and  shape. 

In  this  section  we  will  model  these  perturbations  in  order  to  identify 
initial  tolerances  and  develop  strategies  for  making  the  measurement  in 
an  imperfect  cylinder. 

Let  us  suppose  that  the  deviations  of  the  cylinder  from  perfect  shape 


and  homogeneity  can  be  represented  by  writing  the  cylinder  density  as 


<?(£.<)>)  =  la  sinj<j)+  b  cosjijij  , 

i=0  j=0  L  J  J  J 


[4.44] 


where  bQ  ^  =  1  and  a  is  the  mean  density  of  the  cylinder.  The  P^'s  are 

the  Legendre  polynomials  defined  as  P  =1  and 

o 


P.CC)  =  a2-!)1  . 

1  21i.'  dc1 


[4.45] 


It  is  shown  in  many  standard  texts  on  linear  analysis  that  the  Legendre 
polynomials  form  an  orthogonal  basis  for  the  set  of  piecewise  continuous 
functions  on  the  interval  -1  to  1.  See  for  example  [  Kreider  1966]  .  They 
have  the  additional  property  that  IL  is  even  if  -i  is  even  and  odd  if  -c  is  odd 
Equation  [4.44]  ignores  radial  variations  in  density  since  such  varia¬ 
tions  change  only  in  the  higher  order  terms  in  t/R,  and  these  terms  are  very 
small  as  is  shown  by  [4.12] .  Now  using  [4.13]  and  substituting  the  new 


expression  for  a,  we  have 
„  1  2tt 
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and  substituting  [4.44] 
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Since  we  are  primarily  interested  in  gradients  at  the  center  of  the  cylinder 


we  will  take  z  =  0  and  define  the  error  terms 
1  2tt 
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Since  the  sin(j<j)')  and  cos(j<|>,)  functions  are  a  mutually  orthogonal  set, 
the  only  nonzero  ' s  are  those  with  j  =  0  and  j  =  2.  Also,  since  the 
function  of  4'  in  brackets  is  even,  and  the  is  odd  if  i  is  odd,  only 
E  's  with  even  i  are  nonzero.  Using  these  simplifications 
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are  the  only  nonzero  E..'s. 

ij 

2 

First  we  calculate  some  of  the  E^^g's.  P 2 Cn )  =  %(3  x  -1)  so 


-U  b„  GM 
E20  =  - 2o__ 


(3c,2-i)(2c,2-n2)de' 


[4.51] 
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All  of  these  integrals  can  be  evaluated  explicitly  giving 


to  lowest  order  in  n  .  The  first  eight  nonzero  £^2  anc*  ^2n0  are  listec* 
in  Table  4.1. 

With  these  results  we  can  examine  the  effects  of  various  errors.  The 
error  coefficients  of  the  form  E2n2  are  consi£lerat|ly  larger  than  the  E2n0’ 
but  as  the  gradiometer  is  rotated  in  the  center  plane  of  the  cylinder  the 
gradients  produced  by  the  E errors  have  a  cos2<£  dependence.  Consequently 
if  one  averages  two  gradient  measurements  made  with  the  cylinder  rotated 
by  90°  with  respect  to  the  gradiometer  axis  the  result  is  zero.  Actually 
this  result  is  apparent  from  equation  [4.19] .  If  the  gradiometer  is 
oriented  along  the  z-axis  (the  cylinder  symmetry  axis),  a  mass  distribution 
with  cos2cJ)  dependence  cannot  produce  any  gradient  since  there  is  no  net  mass 
and  all  the  mass  is  equidistant  from  the  gradiometer.  So  from  [4.19]  any 
two  gradients  measured  along  orthogonal  directions  in  a  plane  perpendicular 
to  the  z-axis  must  be  equal  and  opposite.  This  provides  a  simple  way  of 
eliminating  all  errors  of  this  type.  We  simply  measure  horizontal  gradients 
along  the  x  and  y  directions  and  average  them.  Any  difference  in  these 
measurements  would  indicate  the  presence  of  a  $  dependence  in  the  cylinder 
mass  or  density. 

Similarly,  any  of  the  E2ng  errors  can  be  identified  by  comparing  the 


gradient  along  the  z  axis  with  the  sum  of  the  x  and  y  gradients.  In  short, 
if  all  three  gradients  are  measured,  [4.19]  insures  that  their  sum  is  zero 


regardless  of  flaws  in  the  cylinder.  The  cylindrical  geometry  is  still 
crucial,  however,  because  it  relieves  the  experimenter  from  the  need  to 
precisely  position  the  gradiometer  inside  of  the  cylinder. 

4 . 4  Finite  Size  Effects 

Since  the  preceding  calculations  give  the  gradients  at  points,  they 
are  strictly  valid  only  for  inf initesimally  small  gradiometers .  In  this 
section  we  will  consider  the  errors  which  result  from  measuring  gradients 
with  a  finite  size  instrument.  Basically,  we  need  to  know  how  [4.19]  is 
modified  by  the  finite  size  of  the  gradiometer. 

An  indication  that  finite  size  effects  are  small  was  given  in  section 
4.1.  For  example,  we  found  that  the  radial  gradient  changes  by  less  than 
.1%  if  the  gradiometer  is  moved  along  the  z-axis  4.6  cm  from  the  cylinder 
center.  So  if  the  gradiometer  proof  masses  have  a  diameter  less  than  4.6 
cm  we  would  guess  that  the  error  from  its  finite  size  will  be  less  than  .1%. 
The  actual  diameter  of  the  proof  masses  is  4.3  cm. 

Suppose  the  gradiometer  is  oriented  along  the  z-axis  at  the  center  of 
the  cylinder  as  in  4.7a.  To  calculate  the  force  on  the  upper  proof  mass 
we  can  imagine  splitting  the  cylinder  into  two  parts  as  in  4.76.  Since 
the  upper  proof  mass  is  at  the  center  of  the  upper  cylinder  piece,  symmetry 
implies  that  only  the  lower  part  of  the  cylinder  applies  a  net  force  on  this 
part  of  the  gradiometer.  Since  the  gradiometer  is  far  away  compared  to  the 
length  and  thickness  of  this  piece  we  will  calculate  its  gravitational  field 
in  the  vicinity  of  the  gradiometer  using  the  approximation  that  it  is  a  very 
thin  ring  of  material  with  a  mass  per  radian  of  bM/LTT.  Then  the  z-coraponent 
of  acceleration  in  the  center  of  the  upper  cylinder  is 
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[4.59] 


Since  the  gradiometer  is  small  compared  with  R  we  expand  the  denominator  in 
C  =  r/R,  so 
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Integration  over  ((>  gives 
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To  evaluate  the  lowest  order  finite  size  correction  we  integrate  this  accelera 
tion  over  a  disk  of  diameter  d  with  mass  per  unit  area  s,  so 


F  = 
z 


-GMbs 


R 


d/R  2tt 

o 

f  f 

1  15C 

1  I 

2  2l372  f  2  A112 

J  J 

0  0 

^1+r+l/n  2  l+5Z+l/r)1 

[4.62] 


After  carrying  out  the  integration  and  expanding  the  result  in  n  =  2R/L 
one  finds  that 
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where  m  is  the  mass  of  the  disk.  Since  the  gradiometer  is  symmetric  about 
the  origin,  the  result  for  the  lower  disk  is  -F  and  the  gradient  y^  is 


given  by 


gz  =  26m 


1 _  _  3rl"(d/z)"  4  | 

3/2  \  r  21  +  a  ^  '  [ 

I  2|lV]  j 


[4 . 64] 


The  leading  term  is  exactly  the  result  of  section  4.1,  and  the  lowest  order 
correction  is 
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to  lowest  order  in  n  ,  for  a  proof  mass  with  a  radius  of  4.3  cm.  Fortunately 
this  is  quite  a  small  correction  and  does  not  interfere  with  a  .1%  measure¬ 
ment  of  the  cylinder  gradients. 

Next  we  must  consider  the  finite  size  corrections  when  the  gradiometer 
axis  is  perpendicular  to  the  cylinder  axis.  Looking  back  at  equation  [4.5] 
if  we  leave  £  #  0  and  work  out  ar(r,z)  we  find  that  the  first  two  terms  of 
equation  [4.12]  become 


,  ,  1  GM  2  r  T  1  ,  1  if 

(r,z)  =  t-  ri  -  - r  + - r 


From  this  we  can  see  that  the  lowest  order  finite  size  correction  in  q  and 
2 

(2d/L)  comes  from  a  term  of  the  form 


To  calculate  the  gradient  correction  we  need  to  integrate  a'  over  a  disk 
whose  axis  is  along  the  x  axis  and  whose  center  is  a  distance  b  from  the 
origin  and  whose  mass  per  unit  area  is  again  s.  Then 
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This  gives  a  gradient  correction  of 
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This  is  just  half  the  value  of  g'  to  lowest  order  is  n  .  Since  g  /g  =-2, 
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From  this  we  can  see  that  [4.19]  now  does  not  hold  exactly.  Summing  the  finite 
errors  we  have 
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This  effect  puts  a  lower  limit  on  the  value  of  a  which  can  be  detected 
in  equation  [4.36] .  Taking  the  ratio  we  find  that  the  minimum  detectable  a 


is  given  roughly  by 


,  V  *v  _-V 


It  is  interesting  to  consider  briefly  the  finite  size  effect  that 

would  be  produced  by  some  flaw  in  the  cylinder.  Suppose  for  instance, 

that  a  cubic  centimeter  of  cylinder  material  is  missing  in  the  wall 

about  one  meter  from  the  geometrical  center.  It  is  easy  to  calculate 

2 

that  to  the  lowest  order  in  b,  the  baseline,  this  flaw  produces  a  V  <j> 
of  3  x  10  6  E,  an  effect  far  too  small  to  observe  in  this  experiment. 

4 . 5  Experimental  Procedure 

Since  the  d.c.  output  of  the  gradiometer  is  determined  by  the  surround¬ 
ings,  the  signal  from  the  cylindrical  source  must  be  modulated  to  distin¬ 
guish  it  from  the  rest  of  the  gradients  due  to  other  objects.  The  cylin¬ 
drical  geometry  is  convenient  for  this  purpose  since  the  cylinder  can  be 
raised  and  lowered  over  the  cryostat  containing  the  gradiometer. 

Figure  4.8  is  a  rough  sketch  of  the  apparatus  we  will  use  for  moving 
the  cylinder.  In  order  to  support  the  dewar  near  floor  level  we  will  exca¬ 
vate  a  hole  ten  to  twelve  feet  deep  in  the  floor  of  end  station  II  at  the 
high  energy  physics  laboratory.  A  pedestal  will  be  mounted  in  the  center 
of  the  hole  to  support  the  dewar.  The  cylinder  will  be  raised  and  lowered 
with  a  motorized  hoist.  In  the  low  position  the  gradiometer  will  be  at  the 
center  of  the  cylinder,  and  in  the  high  position  the  gradiometer  will  be 
about  10  meters  from  the  center  of  the  cylinder.  To  measure  the  central 
gradients,  the  cylinder  should  be  removed  to  infinity  leaving  everything 
else  fixed.  Since  the  roof  of  the  end  station  prevents  this,  a  correction 
has  to  be  made  for  the  gradient  produced  by  the  cylinder  in  the  high  posi¬ 
tion.  Looking  back  at  [ 4.17]  we  can  calculate  that  the  cylinder  produces 
a  gradient  on  the  high  position  which  is  3.5%  of  the  gradient  in  the  low 
position  if  it  is  moved  two  cylinder  lengths  away  (2z/L  =  4) .  The  signal 
changes  sign  as  the  cylinder  is  moved  away,  so  the  difference  between  the 
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gradient  measured  with  the  cylinder  down  and  the  cylinder  up  has  to  be 
reduced  by  3.5%  to  obtain  the  correct  central  gradient.  This  procedure 
works  for  measuring  the  horizontal  or  vertical  gradient,  and  it  only  needs 
to  be  carried  out  with  an  accuracy  of  a  few  percent  to  obtain  the  central 
gradient  to  an  accuracy  of  .1%. 

By  raising  and  lowering  the  cylinder,  the  signal  frequency  can  be 
moved  away  from  d.c.  and  averaging  can  be  used  to  obtain  a  very  narrow  band 

_Q 

width.  With  the  current  noise  level  of  100E//h¥,  a  bandwidth  of  10  Hz  is 

_2 

necessary  to  obtain  the  desired  noise  level  of  10  E,  so  some  reduction  in 
the  noise  level  is  still  necessary  to  make  this  procedure  practical. 
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CHAPTER  5 


Conclusions 

We  have  developed  a  superconducting  instrument  for  the  measurement 
of  force  gradients.  This  instrument  uses  a  variety  of  properties  of  super¬ 
conductors  to  achieve  very  high  displacement  sensitivity.  The  measurements 
described  in  Chapter  2  demonstrate  that  the  performance  of  the  instrument 
is  very  close  to  the  expected  behavior,  although  the  problem  of  excess  low 
frequency  noise  remains  to  be  solved.  This  problem  is  not  a  consequence  of 
the  techniques  used  in  building  the  instrument,  but  rather  of  the  desire 
to  have  an  instrument  which  can  measure  vertical  and  horizontal  gradients. 

An  instrument  could  be  designed  for  horizontal  use  only  with  higher  compli¬ 
ance  along  the  sensitive  axis.  This  would  reduce  the  gradient  equivalence 
of  a  temperature  change.  However,  as  discussed  in  Chapter  2,  there  are 
steps  which  can  be  taken  to  improve  this  instrument.  The  most  attractive 
of  these  steps  would  be  to  operate  the  instrument  in  a  helium  bath  cooled 
below  the  A  point. 

This  instrument  will  be  incorporated  into  an  experiment  to  check  the 

inverse  square  law  of  gravitation.  The  possibility  of  an  inverse  square  law 

violation  is  briefly  discussed  in  Chapter  3.  In  spite  of  considerable  effort 

by  many  groups  to  verify  the  inverse  square  law,  the  possibility  of  a  large 

3 

violation  with  a  range  between  100  m  and  10  km  or  less  than  1  cm  remains. 

Our  experiment  is  optimized  to  detect  a  small  violation  with  a  range  near 
17  cm,  but  the  use  of  superconducting  technology  for  sensitive  detectors 
may  well  contribute  to  a  variety  of  experiments. 
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